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An  empirically  baaed  mathematical  model  and  computer  code  capable  of  predicting  cumulative  damage  effects  on 
fatigue  crack  growth  in  engine  disks  has  been  developed  for  Waapaloy  (PWA  1007)  and  INI 00  (PWA  1073).  The  com¬ 
puting  software  segregates  a  complex  loading  mission,  applies  a  multi-parameter  interpolation  procedure  to  determine 
cycle-by-cycle  crack  growth,  and  performs  a  specimen  life  calculation.  The  model  was  developed  from  a  broad  matrix 
of  testa  of  both  materials,  and  a  number  of  auxiliary  investigations  were  conducted  to  examine  the  range  of  appli¬ 
cability  of  the  model. 
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SUMMARY 


Cumulative  damage  fracture  mechanics  at  elevated  temperatures  is  particularly  complex, 
requiring  a  basic  understanding  of  the  variables  affecting  crack  growth  and  an  ability  to  model 
material  behavior.  The  synergistic  effects  of  varying  load  and  load  sequencing  on  crack 
propagation  in  engine  disk  materials  preclude  the  use  of  conventional  linear  damage  accumula¬ 
tion  techniques  (e.g.,  Miner’s  Rule)  for  accurate  life  predictions.  Failure  to  consider  interaction 
effects  may  lead  to  significant  errors  in  predicted  crack  growth  behavior  and  in  the  resulting  life 
prediction.  In  most  cases,  the  result  is  an  underestimation  of  propagation  life  and  the  concomi¬ 
tant  economic  penalties  of  premature  removal.  Also  to  be  considered  is  any  overestimate  of  life 
and  the  serious  implications  concerning  component  safety. 

The  objective  of  this  contract  was  to  develop  and  demonstrate  an  empirically  based 
mathematical  model  and  computer  code  capable  of  predicting  cumulative  damage  effects  on 
fatigue  crack  growth  in  gas  turbine  engine  disks. 

Since  the  technical  approach  was  empirical  in  nature,  no  attempt  was  made  to  isolate  and 
describe  individual  microscopic  events  associated  with  propagation  interrupted  by  overloads 
and/or  dwells,  but  rather  synergism  resulting  from  various  combinations  of  these  was  empiri¬ 
cally  described  using  an  interpolative  mathematical  model.  The  resulting  descriptions  of 
generic  elemental  cumulative  damage  events  could  then  be  used  to  analyze  specific  complex 
mission  spectra,  and  the  life  under  such  a  mission  could  then  be  computed  by  simple  numerical 
integration  of  segregated  elements.  A  three-step  process  of  mission  segregation,  crack  growth 
rate  interpolation,  and  calculation  of  fatigue  crack  propagation  life  was  established.  This 
computerized  model  performs  Mission  Analysis  and  Prediction  (MAP)  of  crack  growth  in  gas 
turbine  engine  disks  subject  to  complex  mission  loading. 

Development  of  the  MAP  computing  software  was  based  on  an  analysis  of  representative 
missions  for  TF30  and  F100  engines  employing  Waspaloy  and  INlOO  turbine  disks.  Such 
mission  usage  is  composed  of  a  variety  of  operational  activities  including  aircraft  takeoff,  ferry 
and  refueling,  terrain-following  radar  (TFR),  bombing  runB,  and  combat  under  subsonic  and 
supersonic  conditions. 

Analysis  of  the  turbine  disk  missions  revealed  that  these  low-cycle  fatigue  loading  sequen¬ 
ces  differ  significantly  from  the  high-cycle  loading  spectra  typically  experienced  by  airframe 
components.  The  load  sequence  effects  observed  in  turbine  disk  operation  are  characterized  by  a 
high  occurrence  frequency  of  major  throttle  excursions  (overloads).  As  a  result,  immediate 
post-overload  load  sequence  behavior  is  of  increased  significance,  while  the  long-term  effects  of 
an  overload  on  crack  growth  are  truncated  by  the  frequently  repeated  overload.  The  instan¬ 
taneous  behavior  of  a  crack  propagating  under  fatigue  interrupted  by  frequently  occurring 
overloads  was  shown  to  be  dominated  by  delayed  retardation.  That  is,  both  Waspaloy  and 
INlOO  exhibit  continuous  post-overload  deceleration  in  crack  growth  until  the  recurrence  of  an 
overload  truncates  and  restarts  the  load  interaction  process.  This  behavior  is  not  effectively 
described  by  conventional  models  of  load  sequence  effects.  However,  the  methodology  employed 
in  the  MAP  system  was  developed  specifically  for  prediction  of  mission  crack  growth  with 
frequently  occurring  overloads. 

Development  of  the  Mission  Analysis  and  Prediction  system  was  based  on  the  P&WA 
philosophy  of  empirical  synergistic  modeling.  This  philosophy  states  that  any  complex  mission 
spectrum  can  be  segregated  into  elemental  damage  events  that  can  be  quantitatively  described. 

The  crack  propagation  life  expected  under  such  a  spectrum  can  then  be  computed  as  the  linear 
addition  of  the  damage  associated  with  properly  segregated  events. 

Mission  segregation  is  accomplished  by  the  first  of  three  MAP  algorithms:  the  segregation, 
interpolation,  and  computation  algorithms.  These  are  described  below. 
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The  set  of  rules  by  which  a  mission  profile  is  separated  (segregated)  into  its  elemental 
cumulative  damage  events  constitutes  the  segregation  algorithm.  As  formulated,  the  segrega¬ 
tion  algorithm  translates  an  input  isothermal  stress-time  profile  into  a  parametric,  cycle-by¬ 
cycle  description  of  the  loading  sequence.  This  segregated  mission  definition  provides  input 
compatible  with  the  interpolation  algorithm. 

The  fundamental  strength  of  the  MAP  system  is  the  interpolative  capability  provided  by 
the  hyperbolic  sine  model.  The  application  of  this  procedure  is  executed  by  the  interpolation 
algorithm.  As  formulated,  this  algorithm  performs  a  complex  interpolation  to  obtain  a  crack 
growth  rate  curve  (da/dN  us  AK)  as  a  function  of  the  operating  parameters  (e.g.,  stress  ratio, 
temperature,  cyclic  frequency,  and  load  sequence  parameters).  These  interpolated  curves  are 
subsequently  used  in  the  computation  algorithm  to  obtain  a  fatigue  crack  life  calculation. 

A  cycle-by-cycle  integration  of  the  collection  of  interpolated  crack  growth  curves  is  per¬ 
formed  by  the  computation  algorithm.  The  resulting  life  calculation  is  a  prediction  of  crack 
propagation  resulting  from  the  input  mission  usage  profile. 

The  accuracy  of  the  cumulative  damage  fracture  mechanics  system  was  evlauated  by 
conducting  a  number  of  model  demonstration  tests  simulating  engine  usage.  Using  the  MAP 
system,  specimen  life  calculations  were  obtained  for  the  simulated  missions.  Subsequently, 
fatigue  crack  propagation  tests  were  conducted  to  evaluate  the  model  accuracy. 

The  results  of  the  complete  model  demonstration  program  established  the  utility  of  the 
MAP  system  and  illustrated  the  importance  of  accurate  assessment  of  the  cumulative  damage 
fatigue  crack  growth  which  results  from  mission  usage  of  a  miltary  gas  turbine  engine. 
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SECTION  I 
INTRODUCTION 


Cumulative  damage  fracture  mechanics  at  elevated  temperatures  is  particularly  complex, 
requiring  a  basic  understanding  of  the  variables  affecting  crack  growth  and  an  ability  to  model 
material  behavior.  The  synergistic  effects  of  varying  load  and  load  sequencing  on  crack 
propagation  in  engine  disk  materials  preclude  the  use  of  conventional  linear  damage  accumula¬ 
tion  techniques  (e.g.,  Miner's  Rule)  for  accurate  life  predictions.  Failure  to  consider  interaction 
effects  may  lead  to  significant  errors  in  predicted  crack  growth  behavior  and  in  the  resulting  life 
prediction.  In  most  cases,  the  result  is  an  underestimation  of  propagation  life  and  the  concomi¬ 
tant  economic  penalties  of  premature  removal.  Also  to  be  considered  is  any  overestimate  of  life 
and  the  serious  implications  concerning  component  safety.  The  desirability  of  an  accurate 
crack  propagation  prediction  methodology  is  clearly  evident. 

The  objective  of  this  program  was  to  develop  and  demonstrate  an  empirically  based 
mathematical  model  and  computer  code,  capable  of  predicting  cumulative  damage  effects  on 
fatigue  crack  growth  in  engine  disks. 

Pratt  &  Whitney  Aircraft  Group,  under  the  AFML  sponsorship,  recently  developed  an 
efficient  and  accurate  empirical  method  for  analyzing  synergism  in  crack  propagation  at 
elevated  temperatures  (References  1,  2,  and  3).  No  attempt  was  made  to  isolate  and  describe 
individual  microscopic  events  associated  with  propagation  interrupted  by  overloads  and/or 
dwells,  but  rather  synergism  resulting  from  various  combinations  of  these  was  empirically 
described  using  an  interpolative  mathematical  model.  The  resulting  descriptions  of  generic 
elemental  cumulative  damage  events  were  then  used  to  analyze  specific  complex  mission 
spectra.  The  life  under  such  a  mission  was  then  computed  by  simple  numerical  integration  of 
segregated  elements.  This  proven  method  formed  the  basis  of  our  technical  approach  to  cumula¬ 
tive  damage  fracture  mechanics. 

The  program  consisted  of  a  24-month,  four-phase  technical  effort.  Phase  I,  Mission  Defini¬ 
tion,  defined  representative  composite  missions  for  two  advanced  Air  Force  fighter  engines 
that  employ  Waspaloy  and  INlOO  turbine  disks.  These  results  provided  the  foundation  for  test 
matrix  formulation.  Phase  II,  Mission  Segmentation,  defined,  segregated,  and  quantified  ele¬ 
mental  cumulative  damage  events.  Phase  III,  Mathematical  Model  Development,  evaluated  the 
data  generated  in  Phase  II,  integrated  these  results  with  relevant  pre-existing  data,  and 
translated  the  resulting  model  into  a  computer  code.  Phase  IV  was  Model  Demonstration.  The 
Air  Force  Project  Engineer  provided  two  representative  mission  profiles,  upon  which  life 
history  (a,  N)  predictions  were  made.  These  were  recorded  with  the  Air  Force  Project  Engineer 
prior  to  demonstration  testing.  Fatigue  crack  propagation  (FCP)  tests  on  the  two  materials, 
INlOO  and  Waspaloy,  addressing  the  two  profiles,  and  using  both  compact  and  surface  flaw 
specimens,  were  performed  to  verify  the  model  experimentally.  A  flow  chart  of  the  program 
operating  plan  is  shown  in  Figure  1. 
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Figure  1.  Flow  Chart  of  Program  Operating  Plan 
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SECTION  II 

TECHNICAL  PROGRAM  AND  ACCOMPLISHMENTS 


A.  PHASE  I  •  MISSION  DEFINITION 

1.  Mission  Operating  Profit* 

The  objective  of  the  analysis  of  mission  usage  data,  the  Phase  I  activity,  was  to  provide  a 
basis  for  formulation  of  a  test  and  evaluation  program  for  optimum  model  applicability. 
Mission  usage  field  data  for  military  aircraft  was  generally  obtained  from  pilot  interviews  and 
event  history  recorders  (EHR)  which  have  been  installed  on  advanced  Air  Force  fighters  to 
acquire  real-time  engine  operating  parameters.  The  compilation  and  characterization  of  these 
raw  data  define  representative  engine  usage  for  a  variety  of  operational  activities  including 
aircraft  takeoff,  ferry  and  refueling,  terrain  following  radar  (TFR),  bombing  runs,  and  combat 
under  subsonic  and  supersonic  conditions 

Representative  missions  were  analyzed  for  two  advanced  Air  Force  fighter  engines  employ¬ 
ing  Waspaloy  and  IN  100  turbine  disks.  For  Waspaloy,  a  composite  mission  has  been  prepared 
to  reflect  the  broad  range  of  operating  conditions  imposed  on  this  disk.  This  composite,  pre¬ 
sented  in  Figure  2,  illustrates  the  representative  duty  cycles  in  terms  of  rotor  speed  for  the 
collection  of  individual  mission  usages.  Figure  3  shows  a  similar  mission  for  a  second  engine 
which  employs  turbine  disks  fabricated  from  IN  100.  This  operating  history,  given  in  terms  of 
maximum  stress,  is  based  on  a  Tactical  Air  Command  (TAC)  composite  mission  supplied  by  the 
Air  Force  Materials  Laboratory  in  1974.  This  mission  has  excluded  load  sequences  correspond¬ 
ing  to  operation  of  TFR,  since  this  is  not  a  characteristic  operating  activity  for  this  engine.  The 
two  flight  missions  presented  are  currently  used  as  the  basis  for  field  life  predictions  for 
Waspaloy  and  IN100  disks. 

Cumulative  damage  element  event  definition,  required  specimen  testing,  and  model  devel¬ 
opment  were  based  on  these  two  missions. 

2.  Mission  Stress  Analysis 

A  generalized  fracture  mechanics  residual  life  prediction  system  capable  of  addressing 
crack  propagation  in  disks  under  missiqn  loading  must  consider  all  fracture  critical  locations  in 
the  disk.  These  include:  (1)  bore  (internal  defects),  (2)  rim  (cooling  holes,  attachments),  and 
(3)  web  (boltholes).  Bore  problems  sire  usually  associated  with  inclusions  or  other  intrinsic, 
subsurface  defects  where  crack  propagation  is  not  affected  by  the  surrounding  environment. 
Alternatively,  the  fracture  behavior  of  surface  cracks  in  wrought  nickel-base  superalloys  may 
be  influenced  by  strong  environmental  effects  (e.g.,  oxidation,  sulfidation).  This  program 
addressed  the  propagation  of  surface  cracks  in  air;  the  propagation  of  internal  defects,  or  of 
surface  cracks  subjected  to  extremes  in  environment,  was  not  considered. 

The  stress  analysis  of  fracture  critical  disk  locations  concentrated  on  boltholes  and  cooled 
and  uncooled  rims  of  turbine  disks.  At  these  locations,  stress  histories  depend  on  the  operating 
mission  profile  (Figures  2  and  3),  transient  thermal  gradients,  and  the  geometrical  stress 
concentration.  Under  a  program  for  the  Air  Force  Aero  Propulsion  Laboratory,  “Structural  Life 
Prediction  and  Analysis,”  (Reference  4)  a  survey  was  made  of  estimated  bolthole  stresses  in  the 
two  advanced  Air  Force  fighter  engines  considered  in  the  present  contract.  The  results  of  the 
bolthole  stress  survey  were  examined  as  part  of  the  cumulative  damage  stress  analysis  per¬ 
formed  under  the  current  contract.  In  addition,  bore,  bolthole,  and  rim  stress-temperature-time 
profiles  for  both  Air  Force  and  Navy  turbine  disks  were  reviewed. 


3 


TFR  Bomb  Run 


Figure  2.  Representative  Mission  for  Waspaloy  Disks 


There  are  as  many  different  stress-temperature-time  profiles  as  there  are  fracture  critical 
locations  on  individual  disks.  Some  of  these  are  particularly  complex.  Figure  4  presents  the 
operational  conditions  at  the  fracture  critical  location  for  a  cooled  turbine  rim  (radial  cooling 
hole)  in  a  Waspaloy  disk.  Figure  5  presents  the  stress-temperature-time  response  for  a  bolthole 
in  another  Waspaloy  disk.  The  stress-temperature-time  response  of  a  bolthole  in  an  IN  100  disk 
subjected  to  a  simplified  subsonic  combat  mission  is  given  in  Figure  6. 

The  ultimate  objective  of  this  contract  was  to  produce  results  having  general  applicability 
to  the  analysis  of  surface  cracks  in  Waspaloy  and  IN100  turbine  disks.  Crack  propagation  data 
were  generated  from  laboratory  specimens  having  well-defined  stress  intensity  solutions  and 
thickness  representative  of  disk  operation.  These  data  represent  crack  growth  in  a  uniform 
stress  field;  the  influence  of  aberrations  in  the.elastic  stress  field,  such  as  caused  by  notches,  is 
excluded  from  the  actual  test  data.  The  problem  of  predicting  fatigue  crack  propagation  (FCP) 
in  large  stress  gradients  may  be  solved  analytically  using  advanced  crack  tip  stress  intensity 
solutions  such  as  K  developed  by  P&WA. 


B.  PHASE  II  -  MISSION  SEGMENTATION  AND  TESTING 
1.  Philosophy 

A  thorough  understanding  of  fatigue  crack  propagation  under  complex  cycling  is 
fundamental  to  both  damage-tolerant-design  and  retirement-for-cause  methodologies.  Previous 
investigations  have  considered  certain  aspects  of  FCP  under  streBs-temperature-time  and 
geometry  (thickness)  conditions  representative  of  turbine  airfoil  operation;  however,  additional 
study  of  crack  propagation  under  stress-temperature-time  conditions  characteristic  of  turbine 
disk  operation  was  needed. 

The  goal  of  Phase  II  of  this  contract  was  to  develop  an  understanding  of  cumulative 
damage  FCP  suitable  for  development  of  an  empirical  model  capable  of  accurately  predicting 
crack  growth  for  typical  USAF  turbine  disk  stress-temperature-time  profiles.  The  approach  w  as 
to  employ  a  new  experimental  methodology  which  examines  the  macroscopic  effects  of  realistic 
service  loading  on  FCP.  This  approach,  first  demonstrated  in  AFML  contract  F33615-75-C-5057 
“Applications  of  Fracture  Mechanics  at  Elevated  Temperatures”  (Reference  1,  2,  and  3),  has 
been  extended  and  optimized  for  cumulative  damage  crack  growth. 

The  P&WA  philosophy  of  empirical  synergistic  modeling  is  that  any  complex  mission 
spectrum  can  be  segregated  into  elemental  damage  events  which  can  be  quantitatively 
described.  The  crack  propagation  life  expected  under  such  a  spectrum  can  then  be  computed  as 
the  linear  addition  of  the  damage  associated  with  properly  segregated  events.  The  following 
definition  is  fundamental  to  understanding  this  approach: 

•  An  elemental  damage  event  is  the  smallest  repeating  load-time  sequence 
which  results  in  FCP  not  predictable  by  linear  damage  accumulation 
alone. 

An  analysis  of  the  missions  currently  used  for  field  life  predictions  of  two  advanced  USAF 
fighter  engines  established  cumulative  damage  elemental  events.  These  events  are  defined  as: 
(1)  constant  load  amplitude  cycling,  including  dwells,  (2)  sustained  load  crack  growth,  and 
(3)  overload(s)/LCF  cycling. 


6 


Figure  6.  High-Pressure  Turbine  2nd-Stage  Disk  Bolthole  Stress  vs  Time 


2.  Metallurgical  Analysis 


a.  Waapaloy  (PWA  1007) 

Three  Waspaloy  forgings  were  obtained  from  Reisner  Metals,  Inc.  The  forgings  were 
fabricated  from  Teledyne  Allvac  billet,  heat  No.  L568,  code  WAH.  The  forgings  were  457  mm 
(18  in.)  OD  by  40.6  mm  (1.6  in.)  thick.  The  chemical  composition  of  Waspaloy  (PWA  1007)  met  or 
exceeded  the  specifications  shown  in  Table  1 ,  and  the  results  of  mechanical  property  tests,  also 
given  in  the  table,  satisfied  material  specifications.  Typical  microstructure  of  the  Waspaloy 
forgings  is  shown  in  Figure  7,  and  thin  foil  electron  micrographs  of  representative  structure  are 
presented  in  Figure  8. 

b.  GATORIZED®  IN100  (PWA  1073) 

Two  IN100  forgings  were  GATORIZED  at  P&WA/Florida.  Material  source  was  Udimet 
Powder  Division  of  Special  Metals  Corporation,  heat  codes  F822,  B1  and  B2.  The  mechanical 
properties  and  chemical  composition  of  the  IN100  (PWA  1073)  met  or  exceeded  specifications 
shown  in  Table  2.  Typical  microstructure  of  the  IN100  is  given  in  Figure  9.  Note  the  small  grain 
size,  ASTM  12-14.  Thin-foil  transmission  electron  micrographs  of  representative  dislocation 
substructure  are  given  in  Figure  10. 

3.  Experimental  Program 

a.  Applicability  of  Unaar  Elaatlc  Fracture  Mechanics 

A  fundamental  assumption  implicit  in  the  development  of  the  cumulative  damage  crack 
propagation  model  is  the  applicability  of  linear  elastic  fracture  mechanics  (LEFM)  for 
Waspaloy  and  IN100.  An  examination  of  this  assumption  is  discussed  below.  Consideration  of 
nonlinear  fracture  mechanics  concepts  is  beyond  the  scope  of  this  contract. 

The  presence  of  a  crack  in  a  stressed  component  necessitates  redistribution  of  stresses 
around  the  crack.  The  stress  intensity  factor  is  a  parameter  that  reflects  this  redistribution  and 
is  a  function  of  nominal  stress,  crack  size,  and  specimen  and  crack  geometries.  The  concept  of 
stress  intensity  factor  was  originally  defined  for  an  infinitely  Bharp  elastic  crack  tip,  and  it  is 
this  deformation  that  gives  the  material  resistance  to  crack  propagation. 

The  degree  of  brittleness  of  a  material  (and  the  limit  to  the  applicability  of  linear  elastic 
fracture  mechanics)  is  directly  related  to  the  type  of  stress  redistribution  process  that  occurs  at 
the  crack  tip.  In  the  high  temperature  fatigue  process,  this  redistribution  of  stress  is  expected  to 
depend  on  the  relative  degree  of  elastic,  plastic,  creep,  and  chemical  work  expended  at  the  crack. 
In  a  completely  brittle  material,  relaxation  of  the  crack  tip  stress  field  is  negligible,  and  simple 
reinitiation  of  a  stopped  crack  is  sufficient  to  promote  complete  fracture.  The  absolute  limit  to 
the  applicability  of  fracture  mechanics  is  general  yielding. 

The  usefulness  of  LEFM  depends  on  a  uniparametrical  relationship  between  crack  growth 
rate  and  the  stress  intensity  factor.  Crack  tip  inelasticity,  due  to  material  response  at  elevated 
temperatures,  can  preclude  general  utility  of  K  as  the  correlative  parameter. 

The  following  provides  our  approach  to  the  development  of  a  system  of  empirical  models  for 
accurate  predictions  of  surface  crack  growth  in  Waspaloy  in  IN100  turbine  disks. 

Tests  are  conducted  to  ensure  that  the  crack  tip  experiences  sufficient  geometric  constraint 
and/or  environment  embrittlement  as  to  render  crack  growth  rate  —  AK  relationships 
independent  of  component  geometry  (e.g.,  thickness). 
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TABLE 1 

CHEMICAL  COMPOSITION  AND  MECHANICAL  PROPERTIES  OE  WASPAI.OV 

(PWA  1007) 


Chemical  Composition  (%) 
Minimum  Maximum 


1 


Carbon 

0.02 

0.10 

Manganese 

0.75 

Sulfur 

- 

0.020 

Silicon 

0.75 

Chromium 

18.00 

21.00 

Colbalt 

12.00 

15.00 

Molybdenum 

3.50 

5.00 

Titanium 

2.75 

3.25 

Aluminum 

1.20 

1.60 

Zirconium 

0.02 

0.12 

Boron 

0.003 

0.010 

Iron 

-- 

2.00 

Copper 

- 

0.10 

Bismuth 

- 

0.00005 

Lead 

0.0010 

Nickel 

remainder 

(0.5  ppm) 
(10  ppm) 


Mechanical  Properties 


Tensile 


0.2%  Yield 

Ultimate 

Reduction 

Temperature 

Strength 

Strength 

Elongation 

in  Area 

(  c) 

rn 

(MPa)  (ksi) 

(MPa)  (ksi) 

(%) 

(V 

Room 

Temp 

1062  (154.0) 

1373  (199.2) 

24.5 

37.2 

538 

(1000) 

949  (137.7) 

1251  (181.4) 

17.0 

31.2 

Stress  Rupture 

Time  to 

Temperature  Applied  Stress  Failure  Elongation 
(°C)  (6F)  (MPa)  (ksi)  (hr)  (%) 


732  (1350)  552  (80.0)  56.7 

816  (1500)  310  (45.0)  67.7 


28.5 

27.0 


Mag:  100X 


Mag:  10.000X 


Fn  l  «:119 

Figure  7.  Typical  Microstructure  of  Waspaloy  ( PWA  1007)  {A STM  Grain 
Size  3  to  5) 
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Mag:  10.000X 


Mag:  70.000X 


FP  l:t£t*l 


Figure  8.  Thin  Foil  Transmission  Electron  Micrographs  of  Representative 
Structure  from  a  Waspaloy  (PW’A  1007)  Pancake  Forging 
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I  Alii. I-  L‘ 

CHl.MIC.AI.  COMPOSIl  ION  AN1)  M1.C11AMCA1. 1’ ROl’l.R  l  U.S  CM  C.A  IORI/.I  D® 

IN  1 00  (I’W'A  107.31 


Chemical  Composition  laL) 
Minimum  Maximum 


Carbon 

Manganese 

Sulfur 

Phosphorus 

Silicon 

Chromium 

Cobalt 

Molybdenum 
Titanium 
Aluminum 
Vanadium 
Boron 
Zirconium 
Tungsten 
Iron 
Copper 
Columbium 
and  Tantalum 
Lead* 
Bismuth* 
Oxygen 
Nickel 


•If  determined 

Mechanical  Properties 
'  tensile 


0.2  %  Yield 

Ultimate 

Reduction 

Temperature 
(°C)  (°F) 

Strength  Strength 

(MPa)  (ksi)  (MPa)  (ksi) 

Elongation 

III 

in  A  rea 

Ill 

Room 

Temp 

1166 

(169.1)  1631  (236.5) 

16.5 

29.3 

704 

(1300) 

1113 

(161.4)  1241  (180.0) 

22.0 

27.0 

Stress  Rupture 

Time  to 

Reduction 

Temp 

(°C) 

erature 

Applit 

jMFa) 

'd  S  tress 

Ml 

Failure 

M 

Elongation 

Ill 

in  A  rea 

III 

732 

(1350) 

655 

(95.0) 

22.1 

7.5 

15.6 

732 

(1350) 

655 

(95.0) 

29.0 

15.0 

25.3 

732 

(1350) 

655 

(95.0) 

27.4 

16.7 

24.1 

Creep 

Time  to 

Temp 

TEr 

erature 

Applied  Stress 
iMFal  tfsT) 

Failure 

ML 

Hours  to  Creep 

SS  iSL 

704 

(1300) 

551.6 

(80.0) 

161.0 

127.1 

161.0 

0.05 

0.09 

0.020 

- 

0.010 

0.010 

0.010 

11. 90 

12.90 

18.00 

19.00 

2.80 

3.60 

4.15 

4.50 

4.80 

5.15 

0.58 

0.98 

0.016 

0.024 

0.04 

0.08 

- 

0.05 

- 

0.30 

- 

0.07 

0.04 

0.0002 

(2  ppm) 

0.00005 

(0.5  ppm) 

0.010 

(100  ppm) 

remainder 
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Mag:  10.000X 


FT)  112121 


Figure  9.  Typical  Microstructure  of  GA  TORIZED®  IN  100  (PWA 
(ASTM  Grain  Size  12  to  14) 
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Mag:  70.000X 

Fn  iszia 


Figure  10.  Thin  Foil  Transmission  Electron  Micrographs  of  Representative 
Structure  from  a  GATORIZED®  IN  100  Pancake  Forging 
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While  it  is  known  that  crack  gTowth  rate  is  related  to  stress  intensity  factor  for  specific 
conditions,  tests  are  performed  to  study  the  generality  of  the  K  concept.  In  one  series  of  tests, 
crack  growth  data  are  obtained  in  a  previously  investigated  range,  but  the  specimen  net  section 
stress  is  allowed  to  approach  and  exceed  yield  conditions,  and  deviation  from  elastic  conditions 
is  quantified.  Specimen  geometry  (through  thickness  center  crack  (Figure  11a)  and  compact 
specimen  (Figure  lib))  and  thickness  are  varied  to  test  the  uniqueness  of  da  dN  vs  aK 
relationship  at  a  given  operating  condition. 

Those  geometries  and  thickness  are  used  where  da  dN  vs  aK  is  generally  unique  for  a  given 
operating  condition.  Well-documented  (Reference  5)  specimens  are  used  to  avoid  uncertainties 
concerning  K-calibration.  Thicknesses  are  representative  of  turbine  disks  so  as  to  avoid 
potential  mixed-mode  crack  growth  sometimes  observed  in  thin  sections. 

b.  T Ml  Procedure 

Test  specimens  were  precracked  using  procedures  outlined  in  ASTM  E-399.  Precracking 
was  performed  at  room  temperature  at  a  cyclic  frequency  of  10  or  20  Hz;  anomalous  precracking 
effects  were  easily  recognizable. 

Crack  propagation  testing  was  conducted  on  closed-loop  servocontrolled  equipment 
operated  under  load  control.  Cyclic  tests  were  performed  using  isosceles  triangular  load 
waveforms,  and  specimen  heating  was  provided  by  resistance,  clamshell  furnaces  having 
windows  to  allow  observation  of  crack  growth  at  the  test  temperature. 

Crack  lengths  were  measured  on  both  surfaces  of  the  propagation  specimen  using  a 
traveling  microscope.  This  was  facilitated  by  interrupting  the  cyclic  loading  and  applying  the 
mean  test  load.  This  procedure  held  the  specimen  rigid  while  increasing  crack  tip  visibility.  A 
high  intensity  light  was  used  to  provide  oblique  illumination  to  the  crack  and  further  increase 
crack  visibility.  In  general,  crack  length  measurements  were  taken  at  increments  no  larger  than 
0.50  mm  (0.020  in.).  Crack  length  measurements  are  considered  to  be  accurate  to  within  ±0.025 
mm  (0.001  in.). 

c.  Teat  Program 

(1)  Constant  Load  Amplitude  Cycling 

Crack  propagation  under  constant  amplitude  loading  conditions  is  known  to  be  a  function 
of  the  applied  stress  intensity  range  (within  the  limits  of  linear  elastic  fracture  mechanics).  This 
stress  intensity  range,  AK,  may  be  viewed  as  the  driving  force  for  crack  propagation.  Many 
relationships  have  been  developed  to  correlate  observed  crack  growth  rate  and  the  stress 
intensity  factor.  Paris  and  Erdogan  presented  the  simple  relationship; 

•jfi-  ~  C(AK)“  (1) 

where  C  and  n  are  material  constants.  A  major  drawback  to  equations  similar  to  (1)  is  that  they 
do  not  account  for  the  effects  of  load  sequence.  Additionally,  dependencies  of  crack  growth  on 
temperature  and  loading  frequency  make  the  general  application  of  equations  such  as  (1)  more 
difficult. 
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0.76  cm 
(0.30  in.) 


ill. 430  cm 


a.  Through-Thickness  Center  Crack  Tension  Specimen 


1.27  cm 


b.  Compact  Specimen 

Ft)  II0I96A 

Figure  11.  Test  Specimens  Have  Documented  Fracture  Mechanics  A  nalysis 

18 


A  specific  goal  of  the  constant  load  amplitude  testing  was  to  describe  the  effects  of 
temperature,  stress  ratio  and  loading  frequency  on  crack  propagation  of  Waspaloy  (PWA  1007) 
and  IN10Q(PWA  1073).  The  fundamental  tool  used  to  achieve  this  objective  was  developed  by 
P&WA  under  a  separate  AFML  contract  and  has  been  optimized  under  the  current  program. 
The  mode],  which  provides  an  interpolative  capability  for  crack  growth  rate  as  a  function  of 
several  operating  parameters  (e.g.  temperature,  stress  ratio,  frequency),  will  be  presented  at 
length  in  a  later  section. 

The  following  tasks  were  performed  in  order  to  provide  a  basis  for  an  accurate,  interpolative 
model  for  constant  load  amplitude  FCP;  the  matrix  of  these  basic  propagation  tests  is  given  in 
Table  3.  All  tests  employed  saw  tooth,  constant  load  amplitude  cycling. 

(a)  Heat- to- Heat  Variation 

With  few  exceptions,  specimens  tested  in  this  program  were  taken  from  the  same  heats  of 
material  (Waspaloy  and  IN  100)  and  were  different  from  those  of  existing  data.  The  heats  of 
material  for  this  program  met  P&WA  specifications,  and  tests  were  performed  to  ensure 
compatibility  with  existing  data.  In  order  to  determine  and  quantify  any  significant  variations 
from  heat-to-heat,  test  conditions  were  selected  to  overlap  preexisting  data. 

No  significant  variations  in  crack  propagation  behavior  were  found  among  available  heats 
of  material.  The  extremely  rigid  controls  of  alloy  chemistry  imposed  on  both  IN100  and 
Waspaloy  ensure  consistent  composition.  Processing  and  heat  treatment  specifications, 
supported  by  a  quality  control  inspection  of  microstructure  and  mechanical  properties,  also 
limit  the  performance  variability  of  each  alloy.  As  was  shown  in  Figure  9,  the  characteristic 
microstructure  of  GATORIZED®  IN100  (PWA  1073)  is  extremely  uniform,  having  a  very  small 
grain  size  (ASTM  l2  to  14).  The  crack  growth  behavior  of  this  material  is  correspondingly  well 
behaved.  The  microstructure  of  conventionally  forged  Waspaloy  (PWA  1007),  Figure  8, 
possesses  a  much  larger  characteristic  grain  size  (ASTM  3  to  5),  and  the  allowable  variability  in 
grain  size  is  generally  ASTM  2  to  7.  While  the  heats  of  WAspaloy  evaluated  under  this  contract 
did  not  display  significant  heat-to-heat  variability,  the  potential  for  limited  variation  in  the 
elevated  temperature  cracking  behavior  of  this  alloy  exists  as  a  result  of  gTain  size  flexibility. 

(b)  Effect  of  Loading  Rate  (Frequency) 

Tests  to  determine  the  influence  of  cyclic  frequency  on  elevated  temperature  crack 
propagation  in  Waspaloy  were  performed  under  the  current  program,  and  similar  data  has  been 
generated  for  IN  100  (References  1 , 2,  and  3).  The  range  of  frequencies  examined  (0.00833  Hz  (0.5 
cpm)  to  20  Hz)  was  chosen  so  as  to  reflect  actual  loading  of  a  turbine  disk  as  defined  in  Phase  I. 

(c)  Effect  of  Stress  Ratio 

An  investigation  of  the  effect  of  stress  ratio  (R  =  emmM/amiB)  on  crack  growth  was  conducted 
over  the  range  of  0. 1  <  R  <  0.8  for  a  broad  range  of  temperatures.  The  influence  of  negative  stress 
ratios  was  examined  under  an  auxiliary  investigation  to  be  discussed  later. 

(d)  Effect  of  Temperature 

The  effect  of  operating  temperature  on  subcritical  crack  growth  was  investigated  for 
Waspaloy  under  the  test  program  of  this  contract.  The  effect  of  temperature  for  IN100  had  been 
previously  characterized  (Reference  2).  For  both  materials,  the  range  of  interest  was  427  to 
732°C  (800  to  1360°F),  and  the  environment  was  laboratory  air. 
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A 


(2)  Crack  Propagation  Under  Periodic  Load  Dwell 

This  task  examined  the  crack  propagation  of  W aspaloy  and  IN  100  under  fatigue  cycling  with  a 
dwell  at  maximum  load,  Figure  12.  The  testing,  outlined  in  Table  3,  was  conducted  to  character¬ 
ize  the  effect  of  period  of  load  dwell  and  operating  temperature  on  subcritical  crack  growth. 


Figure  12.  Loading  Waveform  for  Periodic  Load  Dwell  Testing 


(3)  Sustained  Load  Crack  Growth 

The  methodology  used  to  determine  sustained  load  crack  growth  relationship  for  W aspaloy 
and  IN100  was  developed  by  P&WA  under  an  AFML  contract,  and  is  described  in  Reference  2. 
The  test  effort  employed  compact  specimens  and  examined  creep  crack  growth  over  the  range  of 
649  to  732°C  (1200  to  1350°F).  This  work,  discussed  in  detail  in  a  later  section,  resulted  in  an 
interpolative  model  of  sustained  load  crack  growth  as  a  function  of  temperature.  The  test  matrix 
for  this  task  is  included  in  Table  3. 

(4)  Synergistic  Crack  Propagation 

The  majority  of  the  crack  growth  studies  outlined  thus  far  employ  constant  load  amplitude 
fatigue  cycling.  These  tests,  while  necessary,  are  not  completely  representative  of  component 
operating  conditions.  As  discussed  under  Phase  I,  Mission  Definition,  gas  turbine  components 
experience  complex  load  and  temperature  histories.  Crack  growth  relationships  are  compli¬ 
cated  by  load  sequence  effects,  and  the  influence  of  elevated  temperature  operation  on  load 
interaction  phenomena  is  unclear. 
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The  objective  of  the  testing  in  this  task  was  to  establish  possible  synergistic  effects  among 
secondary  parameters  affecting  crack  propagation.  Testing  was  designed  to  investigate  the 
effects  of  low  cycle  fatigue  (LCF)  —  dwell  interactions  at  constant  peak  load,  interactions  of 
major  load  excursions  (overloads)  with  less  severe  cycling  with  variable  stress  ratio,  and 
interactions  of  high  amplitude  sustained  load  on  LCF. 

The  generated  data  provided  a  basis  for  development  of  a  mathematical  model  capable  of 
predicting  cumulative  crack  growth  during  a  complex  stress-temperature-time  spectrum  for 
typical  Waspaloy  and  IN  100  turbine  disks. 

A  typical  operating  mission  consists  of  a  complex  thermal-mechanical  stress  history  such 
as  that  shown  in  Figure  4.  Our  approach  is  to  separate  and  isolate  those  stress-time  histories 
most  suspected  of  having  synergistic  effects  on  crack  growth  in  Waspaloy  and  IN100.  This 
approach  to  model  development  calls  for  isothermal  testing  only,  thus  minimizing  complicated 
and  costly  thermal  mechanical  testing. 

For  each  test  performed  under  the  current  task,  crack  propagation  data  was  obtained  using 
procedures  described  earlier,  and  crack  length  (a)  vs  cumulative  cycles  (N)  data  was  generated. 
See  Appendix  C.  The  actual  test  performed,  defined  in  the  following  tasks,  was  composed  of 
repetitive  segments  containing  an  isolated  load  sequence  (elemental  damage  event). 

(a)  Task  1  —  LCF-Duiell  Interaction  at  Constant  Peak  Stress 

The  purpose  of  Task  1  was  to  investigate  possible  synergistic  effects  of  combined  LCF  and 
dwell.  All  tests  were  performed  with  constant  peak  stress  (load),  and  two  temperatures  were 
investigated,  649  and  732°C  (1200  and  1350°F).  A  schematic  of  the  stress-time  history  is  given  in 
Figure  13.  The  test  matrix  is  presented  in  Table  4. 


Cycles 


(10  cpm) 


ANDwell  = 
td 

R 


Number  of  LCF  Cycles 

N-j  =  10  =  1  min 

Dwell  Time,  Minutes 

N-j  =  20  =  2  min 

pmin/pmax 

N-j  =  40  =  4  min 

R  =  0.1 

FI)  IW42 


Figure  13.  Low  Cycle  Fatigue-Dwell  Interaction  at  Constant  Peak  Stress 
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A  Task  1  mission,  M,  was  defined  as  aNiw|  LCF  cycles  plus  a  2-min  dwell.  As  the  test 
matrix  defines,  ANDjmII  varied  so  that  the  total  time  for  LCF  cyclic  crack  growth  was  less  than, 
equal  to,  and  greater  than  the  dwell  period.  The  testa  investigated  the  interaction  as  a  function 
of  the  ratio  of  LCF  to  sustained  load  crack  growth  within  a  single  mission. 

Crack  growth  relations,  da  dN  vs  .iK,  were  obtained  for  each  of  the  test  conditions  using 
procedures  described  later.  The  relations  defining  crack  growth  will  also  be  discussed  in  a  later 
section. 

(b)  Task  2  —  Variable  Stress  Ratio  Interaction  with  Major  Load  t LCF  —  Overload / 

The  purpose  of  Task  2  was  to  investigate  possible  synergistic  effects  of  major  load  excur¬ 
sions  on  LCF  crack  growth  in  Waspaloy  and  IN  100  The  teat  temperatures  were  649  and  732JC 
(1200  and  1350°F).  A  schematic  of  the  stress-time  history  is  given  in  Figure  14  and  the  test 
matrix  for  Task  2  is  presented  in  Table  4. 


All  Tests  10  cpm 

H  •  n* 


Figure  14.  Variable  Stress  Ratio  With  Major  Load  Excursions 


(c)  Task  3  —  Interaction  of  High  Amplitude  Sustained  Load  with  LCF 

The  purpose  of  Task  3  was  to  investigate  the  possible  synergistic  effects  of  high  amplitude 
sustained  loads  on  LCF  of  Waspaloy  and  IN  100.  However,  initial  construction  of  the  test  matrix 
was  accomplished  without  complete  knowledge  of  the  form  of  typical  mission  profiles  giving 
stress  and  temperature  as  a  function  of  time.  The  missions  (Figures  2  thorugh  6),  defined  under 
Phase  1  of  this  contract,  represent  the  results  of  the  most  current  engine  usage  and  analyses 
conducted  and  describe  conditions  encountered  by  a  disk  operated  in  a  military  gas  turbine 
engine.  A  review  of  these  profiles  shed  question  on  the  ultimate  value  the  overload  dwell-LCF 
tests  proposed  under  this  task.  The  form  of  these  tests  is  illustrated  in  Figure  15. 

The  testing  was  designed  to  measure  the  effect  of  a  dwell  at  overload  upon  subsequent  FCP; 
however,  no  such  overload  dwells  were  observed  in  the  mission  profiles.  Since  all  other  load 
sequences  were  addressed  in  testing  and  modeling,  the  tests  of  this  task  were  replaced  with 
other,  more  pertinent,  tests.  These  include  the  auxiliary  negative  stress  ratio  tests  as  well  as 
selected  tests  designed  to  complement  and  extend  the  original  matrix. 
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td  =  Dwell  Time,  min 
N  =  Number  of  Cycles 

M'  « 

Figure  15.  Interaction  oj  Sustained  High  Load  With  LCF 


(5)  Auxiliary  Investigations 

The  testing  described  in  the  following  sections  was  conducted  to  supplement  the  data 
generated  for  use  in  model  development.  The  results  of  these  investigations  provide  insight  into 
the  nature  of  crack  propagation  under  a  variety  of  test  conditions  which  are  beyond  the  scope  of 
the  model  development  for  this  program.  While  not  actually  used  in  crack  propagation  model¬ 
ing,  the  auxiliary  data  provide  a  means  to  assess  the  limits  of  applicability  of  the  model. 

(a)  Effect  of  Negative  Stress  Ratio 

Fatigue  cycling  which  is  characterized  by  tension-compression  loading  (R  <  O)  may  occur 
in  turbine  disks  under  some  circumstances.  For  example,  thermally  and  mechanically  induced 
strain  gradients  may  produce  local  compressive  stresses  in  a  disk  under  mission  operation, 
Figure  4.  The  negative  R  tests  of  Table  5  were  conducted  in  order  to  evaluate  this  effect  on  crack 
propagation  in  Waspaloy  and  IN100. 

(b)  Effect  of  Net  Section  Stress 

The  through-thickness  center  flaw  crack  propagation  specimen  (shown  in  Figure  11a)  was 
used  to  determine  the  effect  of  Btress  levels  greater  than  0.8  times  the  yield  strength  of  Waspaloy. 
Similar  tests  were  performed  on  IN100  previously.  This  specimen  was  selected  due  to  the  small 
initial  crack  size  (approximately  0.762  mm  (0.030  in.))  obtainable.  The  small  crack  provides  the 
capability  to  obtain  the  necessary  high  net  section  stress  with  corresponding  stress  intensities 
(AK  approximately  55  MPa  ~  (50  ksi  in.))  that  are  comparable  to  existing  linear  elastic  crack 
propagation  curves. 
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(c)  Effect  uf  Prior  Plastic  Deformation 


The  purpose  of  this  task  was  to  investigate  the  influence  of  prior  plastic  cycling  on  subse¬ 
quent  crack  propagation.  Crack  growth  in  material  which  has  previously  been  yielded  is 
representative  of  the  initial  stages  of  propagation  of  a  flaw  which  initiated  in  the  region  of  a 
stress  concentrator.  These  tests  are  outlined  in  Table  5. 

(d)  Effect  of  Overload ■  Underload-LCF  Sequencing 

The  crack  propagation  testing  upon  which  the  cumulative  damage  modeling  was  based  is 
isothermal  in  nature.  It  was  assumed  that  a  nonisotherm al  mission  may  be  suitably  repre¬ 
sented  by  a  series  of  isothermal  segments.  The  purpose  of  this  investigation  was  a  brief 
examination  of  crack  propagation  of  Waspaloy  (PWA  1007)  subject  to  thermal-mechanical 
fatigue  (Table  5). 

(e)  Thermal-Mechanical  Fatigue 

The  crack  propagation  testing  upon  which  the  cumulative  damage  modeling  was  based  is 
isothermal  in  nature.  It  was  assumed  that  a  nonisothermal  mission  may  be  suitably  repre¬ 
sented  by  a  series  of  isothermal  segments.  The  purpose  of  this  investigation  was  a  brief 
examination  of  crack  propagation  of  Waspaloy  (PWA  1007)  subject  to  thermal-mechanical 
fatigue  (Table  5). 

(f)  Effect  of  Thickness 

The  goal  of  this  task  was  to  investigate  the  effect  of  specimen  thickness  on  crack  propaga¬ 
tion  in  Waspaloy  (PWA  1007).  Of  specific  concern  was  the  qualification  of  crack  propagation 
data  as  thickness  independent  (i.e.,  representative  of  the  most  severe  geometric  constraint). 


C.  PHASE  III  -  MATHEMATICAL  MODEL  DEVELOPMENT 
1.  Model  Description 

To  analyze  a  complex  mission  spectrum  in  terms  of  consistently  accurate  cumulative 
damage  accounting  requires  two  basic  tools:  (1)  a  mathematical  model  describing  crack 
propagation  at  frequencies  ,  temperatures,  and  stress  ratios  representative  of  turbine  disk 
operation,  and  (2)  a  model  which  describes  cumulative  damage  synergistic  effects  such  as: 

•  Overloads  of  varying  degree  and  occurrence  frequency 

•  Dwells  of  varying  duration 

•  Combinations  of  these  with  simple  cyclic  crack  growth  behavior. 

The  requisite  tool  is  the  Interpolate  Synergistic  Hyperbolic  Sine  Model,  SINH. 


) 
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a.  Mathematical  Formulation 


The  model  is  based  on  the  hyperbolic  sine  equation, 

log  (da  dN>  =  C,  sinh  (C2  (log  UK)  ♦  C,))  +  C4  (2) 

wher  ■  .  le  coefficients  have  been  shown  (References  1,  2,  and  3i  to  be  functions  of  test  fre¬ 
quency,  stress  ratio,  and  temperature: 

C,  -  materia]  constant 

C,  =  f,  (v,  R,  T) 

C‘  =  f*  (f,  R,  T) 

C4  =  f4  (v,  R,  T) 

Additionally,  Equation  2  may  also  be  used  to  describe  crack  propagation  data  as  a  func¬ 
tion  of  load  sequence  parameters. 

Descriptions  of  the  mathematical  model  (SINH)  and  the  hyperbolic  sine  function  (sinh) 
on  which  it  is  based  are  provided  in  the  following  paragraphs. 

b.  Charadarlttlc*  of  the  tlnh 

The  hyperbolic  sine  is  defined  as 

y  =  sinh  x  = — -  (3) 

and  when  presented  on  Cartesian  coordinates  it  appears  as  shown  in  Figure  16.  The  function 
is  zero  at  x  =  0  and  has  its  inflection  there. 

The  introduction  of  the  four  regression  coefficients,  Cj  through  C4,  permits  relocation  of 
the  point  of  inflection  and  scaling  of  both  axes.  In  the  equation 

(y  -  C4)  =  sinh  (x  +  C3),  (4) 

Cj  establishes  the  horizontal  location  of  the  hyperbolic  sine  point  of  inflection  and  C4  locates 
its  vertical  position. 

To  scale  the  axes,  C,  and  C2  are  introduced 

7T~*~ —  =  sinh  (C2(x  +  CJ)  <5> 

which  can  be  rewritten  as 

y  =  C;  sinh  (C2  (x  +  Cj))  +  C4 

of  which  Equation  2  is  a  special  case  where  y  =  log  (da/dN)  and  x  =  log  UK).  Note  that  C3 
has  units  of  log  UK)  and  C4  has  units  of  log  (da/dN);  C,  and  C2  are  dimensionless  and  can 
be  conceptualized  as  stretching  the  curve  vertically  and  horizontally,  respectively.  Expe¬ 
rience  indicates  that,  for  a  given  material,  Ct  can  be  fixed  without  adversely  affecting  model 
flexibility. 
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V 


figure  16.  Hyperbolic  Sine  on  Cartesian  Coordinates 


Experience  has  also  shown  that  the  materials  considered  here  (Waspaloy  and  IN100) 
exhibit  symmetrical  da/dN  vs  aK  relationships  for  simple  cycling.  For  materials  which  do 
not  display  this  symmetry,  an  asymmetric  variation  of  the  model  is  available.  Here,  the 
shape  coefficient,  C2,  is  separated  into  similar  coefficients  C2  and  C'2.  The  first  of  these  coef¬ 
ficients  determines  model  curvature  above  the  point  of  inflection,  while  C'2  controls  curva¬ 
ture  below  inflection.  The  model  remains  a  continuous  function,  and  when  the  data  is  sym¬ 
metric  C2  =  C'T  and  the  model  reverts  to  its  original  form. 
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The  hyperbolic  sine  equation  was  selected  as  the  model  for  the  following  reasons: 

•  It  exhibits  the  overall  shape  of  typical  da  dN  vs  plots,  obtained  over 
several  decates  of  crack  growth  rates. 

•  All  or  part  of  the  equation  may  be  used  to  fit  data,  since  the  sinh  has 
both  a  concave  and  a  convex  half  and  a  nearly  linear  portion  near 
inflection.  Also,  the  slope  at  inflection  can  vary  with  the  fitting  con¬ 
stants.  (By  comparison,  the  slope  of  an  x 1  model  is  always  zero  at  in¬ 
flection.) 

•  The  sinh  is  not  periodic  (e.g.,  trigonometric  tangent)  nor  asymptotic 
(e.g.,  tangent,  or  inverse  hyperbolic  tangent);  therefore,  when  extrapola¬ 
tion  becomes  necessary,  the  sinh  behaves  well  at  distances  removed 
from  the  data,  quite  unlike  most  polynomial,  periodic,  or  asymptotic 
functions. 

•  This  model  requires  no  information  other  than  (a  vs  N)  data.  By  com¬ 
parison,  some  other  models  in  current  use  require  both  Kth  and  Kk,  in 
addition  to  (a  vs  N)  data,  to  model  crack  growth  behavior.  Both  Kth  and 
K[c  are  difficult  to  obtain  experimentally,  Kth  because  of  the  extremely 
small  crack  growth  measurements  necessary  and  KIc  because  of  gross 
plasticity  at  the  crack  tip  encountered  in  fracture-toughness  testing  at 
elevated  temperatures. 

c.  SINH  Description s  ol  Beale  Propagation 

The  hyperbolic  sine  model  is  easily  adapted  to  describe  fundamental  parametric  effects  of 
operating  conditions  such  as  stress  ratio,  frequency,  and  temperature  on  crack  growth  rate. 

Figure  17  schematically  depicts  the  qualitative  effects  on  crack  propagation  rate  of  fre¬ 
quency  (Figure  17a),  stress  ratio  (Figure  17b),  and  temperature  (Figure  17c). 

Because  of  the  simple  relationships  observed  between  the  coefficients  of  the  SINH  model 
and  the  fundamental  propagation-controlling  parameters,  interpolations  are  straightforward. 
It  is  here  that  the  model  demonstrates  its  great  usefulness:  the  hyperbolic  sine  model  provides 
descriptions  of  crack  propagation  characteristics  where  data  are  unavailable.  The  interpolation 
algorithm  is  described  under  Algorithm  Definition  and  Development. 

2.  Advanced  Regression  Considerations 

Interpolative  modeling  of  crack  propagation  as  a  function  of  operating  parameters  such 
as  frequency,  stress  ratio,  and  temperature  (Figure  17)  requires  a  multiple  regression  capabil¬ 
ity  which  allows  simultaneous  consideration  of  several  different  collections  of  data,  each 
differing  from  the  others  by  only  one  FCP  —  controlling  parameter.  Separate  regression  of 
the  individual  data  sets  contributing  to  an  interpolative  model  does  not  allow  data  at  one 
condition  to  influence  the  model  at  another  condition.  However,  the  final  model  was  to  have 
broad  interpolative  capability  with  behavior  at  one  condition  used  to  describe  FCP  at 
another  condition.  Therefore,  it  is  desirable  to  permit  the  data  to  exhibit  their  mutual  influ¬ 
ence  during  the  modeling  process. 


30 


Pratt  &  Whitney  Aircraft  has  developed  a  mathematical  technique  to  accomplish  this. 
Individual  sets  of  data  are  treated  independently  relative  to  some  of  the  SINK  coefficients, 
while  the  entire  collection  is  treated  as  an  entity  with  respect  to  the  interpolative  coefficients 
(functions  of  v,  R,  T,  etc.).  This  improved  method  was  used  to  develop  the  propagation  models 
presented  in  this  report.  The  actual  computational  procedure  employed  in  the  modeling, 
referred  to  as  the  Method  of  Least  Squares,  is  described  in  the  following  paragraphs. 

The  goal  of  this  procedure  is  to  determine  model  coefficients  for  which  the  resulting 
curve  through  the  data  will  have  the  least  (minimum)  summed  squared  error  between  calcu¬ 
lated  and  observed  values  for  the  dependent  variable  (see  Figure  18).  In  this  instance,  the 
independent  and  dependent  variables,  x  and  y,  are  log  (AK)  and  log  (da/dN),  respectively. 


FD  130630 

Figure  18.  Method  of  Least  Squares 
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Define  the  sum  of  the  squared  errors  as 


n 

=  2  E?  =  £  <y«,,  •  y,)2 

i  =  1 


(7) 


Since  yc„  =  f  (C2,  C3,  C4,  x^,  E  is  also  a  function  of  C2,  C3,  C4. 

Now,  E2  will  be  a  minimum  when  each  of  its  partial  derivatives  is  zero  simultaneously. 
That  is 


3E2 

2E3E 

=  0 

(8) 

lc7 = 

ac2 

3E2 

2E3E 

=  0 

(9) 

ac3  - 

ac3 

aE2 

2E3E 

=  0 

(10) 

dCd  - 

dCd 

when  f  is  the  SINH  model, 

Ei  =  C1SINH(C2(xi  +  C3))  +  C4-yi  (11) 

and 


aE 

-gjp-  =  C,  cosh  (C2(xi  +  C3))(xj  +  C3) 
ux-'2 

(12) 

dE 

^  -  C,  cosh  (C2(x,  +  CjJKCj) 

(13) 

aE 

ac«  _1 

(14) 

Now,  substituting  Equations  11, 12, 13,  and  14  into  Equations  8,  9,  and  10,  and  solving*  the 
resulting  three  simultaneous  equations  provides  the  values  for  C2,  C3,  and  C4  for  which 
Equation  7  will  be  a  minimum. 


The  foregoing  discussion  explains  the  procedure  for  determining  the  coefficients  for 
SINH  curve.  Suppose  further  that  each  SINH  representation  of  FCP  is  related  to  the  others, 
the  relationship  depending  on  differences  in  frequency,  stress  ratio,  temperature,  etc.  Con¬ 
sider  frequency  as  an  example,  and  assume  that  the  points  of  inflection  are  linearly  related, 
viz: 


Cj,  j  s  Cjj  ♦  C*  <C4j) 
for  j  different  SINH  curves. 


(15) 


*Io  this  butane*,  the  resulting  ainmltaneou*  aquation*  ar*  nonlinaar  In  Cr  C9  and  C4.  Th*  *olation  therefor* 
r*quin*  *om*  itarati  a  procedure  inch  as  an  N-dl mansions!  Newton- Raphaon  method. 
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Assume  also  that  C2  and  C4  are  related  to  test  frequency  by  Equations  16  and  17: 


C„,  1  =  C;«  +  C.«5  <lo&  M) 

(16) 

C  >.  j  =  C37  +  (log  (l)) 

(17) 

Coefficients  C33  through  can  be  then  determined  by  substituting  Equations  15,  16,  and  17 
into  Equation  11  and  differentiating  with  respect  to  CM  through  in  a  manner  analogous 
to  that  used  in  determining  C2,  C3,  C4  in  the  foregoing  discussion. 

3.  Algorithm  Definition  and  Development 

The  computer  code  developed  under  Phase  III  has  as  its  mathematical  formulation  three 
algorithms:  the  segregation,  interpolation,  and  computation  algorithms. 

a.  The  Segregation  Algorithm 

The  set  of  rules  by  which  a  mission  profile  is  separated  (segregated)  into  its  elemental 
cumulative  damage  events  constitutes  the  segregation  algorithm.  As  formulated,  the  segre¬ 
gation  algorithm  translates  an  input  isothermal  stress-time  profile,  such  as  illustrated  in 
Figure  5,  into  a  parametric,  cycle-by-cycle  description  of  the  loading  sequence.  This  segre¬ 
gated  mission  definition,  given  in  terms  of  the  elemental  mission  segments  addressed  in  the 
test  program,  provides  input  compatible  with  the  interpolation  algorithm. 

Development  of  the  segregation  algorithm  was  based  on  the  mission  definition  described 
under  Phase  I,  Figures  2  through  6.  A  characteristic  of  these  representative  operational  pro¬ 
files  is  a  high  occurrence  frequency  for  major  load  excursions  (overloads).  As  discussed  in 
Appendix  A  (Reference  6),  the  implication  of  this  observation  is  to  diminish  the  importance 
of  long  term  overload  effects  and  to  place  increased  significance  on  the  synergistic  interac¬ 
tion  which  occurs  immediately  following  the  overload.  As  a  result,  a  primary  investigation 
conducted  under  the  test  program  was  to  characterize  the  short  term  effects  of  overloading. 
Similarly,  the  interpolation  algorithm,  developed  to  address  mission  profiles  representative 
of  turbine  disk  operation,  is  not  routinely  applicable  to  missions  which  may  contain  infre¬ 
quent  overloads  (e.g.  airframe  loading  spectra). 

The  general  procedure  employed  in  the  segregation  algorithm  is  outlined  as  follows: 

1.  Define  isothermal  operating  temperature 

2.  Locate  load  dwells 

3.  Locate  major  load  excursions  (overloads) 

4.  Calculate  the  number  of  cycles  between  successive  overloads 

5.  Calculate  overload  ratio  for  each  less  severe  cycle  following  an  overload 

6.  Determine  stress  ratio  and  frequency  of  all  fatigue  cycles. 

Upon  completion  of  mission  segregation,  the  data  is  passed  to  the  interpolation  algorithm, 
and  the  appropriate  crack  growth  rate  curves  are  obtained. 

b.  The  Interpolation  Algorithm 

The  fundamental  strength  of  the  hyperbolic  sine  model  is  its  interpolative  capacity.  The 
procedure  known  as  the  interpolation  algorithm  for  calculating  the  SINH  coefficients,  describ¬ 
ing  FCP  under  the  influence  of  a  segregated  elemental  event,  is  illustrated  in  the  following 
paragraphs.  For  simplicity,  the  example  considers  only  simple  cycling  with  varying  frequency, 
stress  ratio,  and  temperature  (v,  R,  T).  The  actual  computer  software  includes  other  life¬ 
controlling  parameters:  period  of  load  dwell,  overload  ratio,  and  cycles  between  overloads. 
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The  coefficient  (e.g.,  C2,  C3,  and  Ct)  at  any  intermediate  value  of  an  element  life-controlling 
parameter,  can  be  determined  from  Equation  18. 

j=  2,3,4  (18) 


=  interpolated  values  of  coefficients 


=  differences  from  baseline  values. 


Since  the  SINH  coefficients  are  linear  functions  of  the  controlling  parameters,*  it  is  evident 
that: 

[AC2  “I  I"  dC2/dv,dC2/dR,dC2/&T 
AC3  1=1  ac3/dv,  dC/dR,  dC/dT 
AC4  J  L  dC/dv,  dC/dR,  dCt/ar 

NX1  NXN 

where: 

Au  “| 

AR  I  =  differences  from  baseline  values 
AT  J 

and  the  NXN  partial  derivative  matrix  is  easily  determined  from  the  slopes  of  the  lines  relat¬ 
ing  each  coefficient  with  each  rate  controlling  parameter.  The  computation  of  the  intermediate 
coefficients,  using  Equation  18,  is  then  straightforward. 

c.  Tha  Computational  Algorithm 

The  product  of  the  interpolation  algorithm  is  a  collection  of  SINH  descriptions  of  crack 
propagation  corresponding  to  the  segregated  mission  cycles.  The  interpolated  crack  growth 
curves  are  representative  of  the  cycle-by-cycle  crack  growth  under  the  individual  loading  condi¬ 
tions,  and  the  effects  of  load  sequence  synergism  are  implicit.  The  method  of  computation  of 
crack  growth  under  the  complex  mission. 

A  simple  cycle-by-cycle  integration  is  used  to  sum  the  incremental  crack  advances,  Aa^ 
which  comprise  the  incremental  cyclic  life,  AN,,  (or  event  life,  AEj) 


b 

*  Ad  1 

A  R  I 

l 

.  AT  J 

N  X  1 


da/dN  =  f(a,  AK. . .) 

(20) 

Aa,  =  f(a,  AK  . . .)  AN. 

(21) 

•Strictly  ipeaking,  the  coefficient*  are  nonlinear  function*  of  v,  R,  and  T;  however,  they  are  linear  function*  of  other 
function*.  Thi*  aimplification  wa*  used  here  for  presentation  clarity. 
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The  sum  of  all  these  incremental  crack  advances  results  in  the  crack  propagation  during  one 
mission,  ia.^ 


Aa_ 


L 

i=i 


A  a, 


where  n  is  the  number  of  elemental  events  comprising  the  mission  profile. 


The  total  mission  lifetime,  NmiMiOM,  is  similarly  calculated  as  the  number  of  mission  crack 
advances,  Aa  required  to  reach  critical  cracklength. 


N_ 


roiMion* 


Aa_ 


i=l 


This  summed  numerical  integration  procedure  is  referred  to  as  the  computational 
algorithm. 


4.  SINH  Descriptions  ot  Crack  Propagation 

Interpolative  crack  propagation  models  developed  from  the  data  generated  under  the  test 
programs  outlined  in  Tables  3  and  4  are  presented  in  the  following  pages.  Each  of  the  indi¬ 
vidual  models  describes  the  influence  of  a  single  test  variable  (e.g.,  temperature,  frequency, 
stress  ratio)  upon  crack  propagation.  The  combination  of  these  models  forms  a  unified  des¬ 
cription  of  the  crack  growth  behavior  over  the  range  of  operating  conditions  addressed  in  the 
test  matrices.  This  unified  description  is  formalized  in  the  interpolation  algorithm. 

a.  Watpaloy  (PWA  1007)  Crack  Propagation 

The  majority  of  the  Waspaloy  data  presented  in  this  report  were  generated  under  the 
current  contract  and  represent  crack  propagation  in  a  single  heat  of  material.  A  limited 
amount  of  additional  data  (Reference  7),  generated  from  a  second  heat  of  this  same  alloy, 
was  also  used.  The  inclusion  of  additional  data  aided  in  model  refinement  and  extension. 

Regression  of  the  crack  length  vs  cycle  (a,  N)  data  to  produce  crack  growth  rate  vs 
applied  stress  intensity  range  (da/dN,  AK)  data,  in  general,  was  accomplished  with  the 
seven-point  incremental  polynomial  technique  (Reference  8).  This  reduction  technique  was 
found  desirable  for  Waspaloy  data  but  less  beneficial  when  reducing  IN100  data,  since  typi¬ 
cal  scatter  in  (da/dN,  AK)  data  in  Waspaloy  is  considerably  more  severe  than  in  IN100. 

The  severity  of  data  scatter  in  these  two  alloys  is  closely  linked  to  grain  size.  Elevated 
temperature  fatigue  of  Waspaloy  and  IN100  results  in  discontinuous  advances  of  the 
observed  surface  crack,  and  the  increments  of  discontinuous  advance  are  on  the  order  of  the 
material  grain  size.  As  the  magnitude  of  the  discontinuous  crack  advance  (Aa^)  increases, 
the  variability  of  the  increment  in  crack  growth  (Aa)  increases  by  a  value  within  the  range 
±2  (Aadi>c),  and  the  associated  scatter  in  da/dN  increases.  The  large  grain  size  of  Waspaloy 
(ASTM  2-7)  makes  this  a  significant  effect,  while  the  fine  crystalline  structure  of  IN  100 
(ASTM  12-14)  results  in  much  less  inherent  scatter  in  crack  growth  rate. 
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The  mathematical  procedure  employed  in  the  seven-point  incremental  polynomial  tech¬ 
nique  for  data  reduction  produces  five  fewer  (da/dN,  AK)  data  points  than  does  the  direct 
secant  method.  Generally,  this  presents  no  significant  problem;  however,  (a,  N>  data  sets 
containing  a  small  number  of  points  may  be  severely  diminished  during  data  reduction. 
Therefore,  small  data  sets  were  reduced  by  the  direct  secant  method  and  their  specimen 
numbers  bear  the  prefix  “100.”  Data  reduced  with  the  seven-point  incremental  polynomial 
method  have  identification  numbers  beginning  with  “7AN”  (See  Appendix  C). 

(1)  Basic  Propagation,  Constant  Load  Amplitude  Cycling 

(a)  Effect  of  Loading  Rate  (Frequency) 

Experience  with  turbine  disk  alloys  indicates  that  changing  test  frequency,  while  hold¬ 
ing  stress  ratio  and  temperature  constant,  produces  crack  growth  curves  similar  in  shape  but 
shifted  along  a  nearly  vertical  line  passing  through  the  points  of  inflection.  The  location  of 
these  inflection  points  is  related  to  test  frequency;  reduced  frequency  (increased  cycle  dura¬ 
tion)  is  observed  to  e  ccelerate  crack  growth  rate. 

Figures  19  and  20  illustrate  the  effects  of  cyclic  frequency  on  FCP  in  Waspaloy,  at  649°C 
(1200  °F),  R  =  0.1.  The  four  curves  represent  crack  growth  under  frequencies  of  0.00833,  0.167, 
0.333,  and  20  Hz  (0.5  cpm,  10  cpm,  20  cpm,  20  Hz).  Note  that  points  of  inflection  describe  a 
straight  line,  i.e.:  C3  =  CM  +  x  C4.  Figure  21  demonstrates  the  interrelationships  of  the 
other  SINH  coefficients  and  test  frequency.  Note  that  the  coefficients  exactly  describe  a 
straight  line,  so  that  any  aberrations  in  the  model  appear  when  the  SINH  curve  is  plotted 
with  the  data  it  represents.  This  model  completely  describes  the  effect  of  frequency  on  FCP 
in  Waspaloy  at  this  temperature  by  interpolation  of  the  equations  given  in  Figure  21. 

(b)  Effect  of  Stress  Ratio,  0.167  Hz  (10  cpm) 

The  influence  of  stress  ratio  (R)  on  crack  propagation  in  Waspaloy  (PWA  1007)  fatigued 
at  0.167  Hz  (10  cpm)  649°C  (1200°F)  is  illustrated  in  Figures  22  through  24.  Data  from  tests 
at  four  difference  stress  ratios  were  used  to  develop  a  model  (Figures  22  and  23)  of  this  effect 
for  0.05  <  R  <  0.80.  A  plot  of  the  correlative  parameters  and  the  defining  equations  is  given 
in  Figure  24.  Each  of  the  SINH  coefficients  is  defined  as  a  linear  function  of  log  (1  -  R)  for  the 
range  of  positive  stress  ratios. 

This  relationship  should  not  be  extrapolated  into  the  range  of  negative  stress  ratios  (i.e., 
tension-compression  fatigue).  Data  gathered  from  tests  at  R  =  -0.5  and  R  =  -1.0  demonstrate 
the  invalid  nature  of  this  extrapolation.  Crack  propagation  under  negative  stress  ratio 
fatigue  is  discussed  under  “Auxiliary  Investigations”  later  in  this  report. 

(c)  Effect  of  Stress  Ratio,  20  Hz 

The  effect  of  stress  ratio  on  fatigue  crack  propagation  of  Waspaloy  is  a  function  of  the 
associated  cyclic  frequency  of  the  loading.  In  order  to  address  this  functional  dependence, 
testing  was  conducted  at  a  frequency  of  20  Hz  and  0.05  <  R  <  0.80.  The  resulting  model  is 
illustrated  in  Figures  25  and  26.  A  plot  of  the  correlative  parameters  is  shown  in  Figure  27. 
The  combination  of  the  20  Hz  stress  ratio  model  and  the  similar  0.167  Hz  model  provides  full 
interpolative  capability  for  crack  growth  rate  as  a  function  of  stress  ratio  and  frequency 
simultaneously. 
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Figure  19.  Waspaloy  (PWA  1007)  Crack  Propagation,  Frequency  Model, 
649°C(1200°F) 
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Figure  23.  Statistics  for  Waspaloy  (PWA  1007)  Crack  Propagation,  Stress  Ratio  ! 
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Figure  26.  Statistics  for  Waspaloy  (PWA  1007)  Crack  Propagation,  Stress  Ratio  Model,  20.0  //:,  OC’C  ( l  JOO ” !■  j 
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Figure  27.  Waspaloy  ( PWA  1007)  Stress  Ratio  Model  Correlative  Para¬ 
meters,  20  Hz,  649°  C  (1200° F) 
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(d)  Effect  of  Temperature 

The  effect  of  operating  temperature  on  crack  propagation  of  Waspaloy  (PWA  1007)  is 
illustrated  in  the  composite  plot  of  Figure  28.  (See  also  Figure  29.)  This  SINH  model  des¬ 
cribes  the  detrimental  influence  of  increasing  temperature  on  crack  growth  in  this  alloy.  The 
relationships  defining  the  associated  SINH  coefficients  are  presented  in  Figure  30. 

(2)  Crack  Propagation  Under  Periodic  Load  Dwell 

(a)  Effect  of  Load  Dwell,  649° C  (1200°F) 

Crack  propagation  in  Waspaloy  (PWA  1007)  under  periodic  load  dwell.  Figure  12,  at 
649°C  (1200°F)  is  a  function  of  the  period  of  dwell.  Tests  conducted  with  dwells  ranging  from 
120  to  900  sec,  were  used  to  develop  the  interpolative  SINH  model  of  the  effect  of  dwell 
shown  in  Figure  31.  (See  also  Figure  32.)  While  experience  has  demonstrated  that  the  SINH 
coefficients  are  linear  functions  of  the  log  of  the  period  of  dwell,  this  relationship  does  not 
completely  describe  the  behavior  of  data  of  Figure  31.  Coefficients  C2  and  C3  are  linear  func¬ 
tions,  but  two  linear  segments  are  required  to  define  the  behavior  of  C4  with  period  of  load 
dwell.  This  is  shown  in  the  plot  of  Figure  33. 

(b)  Effect  of  Load  Dwell,  732° C  (1350°F) 

The  crack  propagation  behavior  of  Waspaloy  (PWA  1007)  under  periodic  load  dwell  dis¬ 
plays  a  significant  dependence  on  operating  temperature.  A  SINH  model  developed  from 
data  generated  in  Waspaloy  at  732°C  (1350°F)  under  dwells  of  120  and  900  sec,  is  presented 
in  Figure  34.  (See  also  Figure  35.)  The  interpolative  functions  describing  this  behavior  are 
given  in  Figure  36. 

Comparison  of  the  649°C  dwell  data  with  data  generated  under  similar  loading  at  732°C 
indicates  that  crack  propagation  rate  under  900  sec  dwell  conditions  is  nearly  independent  of 
temperature.  However,  under  120  sec  dwell  loading,  crack  propagation  at  732°C  is  signifi¬ 
cantly  slower  than  at  649°C  as  shown  in  Figure  37.  The  reduced  crack  growth  rate  at  the 
higher  temperature  is  attributed  to  crack  tip  blunting  which  results  from  an  appreciable 
creep  component  at  732°C.  During  the  732°C,  900  sec  dwell  test,  crack  tip  blunting  resulted 
in  crack  arrest  at  values  of  AK  <  35  MPa  \fm  and  in  general,  severely  complicated  data 
generation  at  this  condition. 

The  interactive  effect  of  periodic  load  dwell  and  operating  temperature  is  included  in  the 
general  interpolation  algorithm  developed  under  this  contract. 

(3)  Sustained  Load  Crack  Growth 

Conventional  experimental  methods  for  determination  of  crack  propagation  behavior 
under  sustained  loading  call  for  incremental  crack  length  measurements  over  the  duration  of 
the  test.  This  allows  generation  of  crack  growth  rate  vs  applied  stress  intensity  (da/dt  vs  K) 
data  which  characterize  crack  growth  under  the  conditions  of  testing. 

We  have  applied  an  alternative  procedure  to  elevated  temperature  sustained  load  crack 
propagation  testing  (Reference  2).  The  underlying  assumption  of  this  method  is  that  the 
sustained  load  crack  growth  curve  exhibits  a  sigmoidal  shape  definable  by  an  equation  of 
the  form: 

da/dt  =  10  **  (C,  SINH  (C2  (log  K  +  Cj))  +  C4).  (24) 
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Figure  33.  Waspaloy  (PWA  1 007)  Dwell  Model  Correlative  Parameters,  R  - 
0.05,  649° C  (1200° F) 
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COEFFICIENTS  C2  AND  C4  VS.  DWELL  TiMt 
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Figure  36. 


Waspaloy  (PWA  1007)  Dwell  Model  Correlative  Parameters, 
R  =  0.10.  732° C  (1350° F) 
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The  experimental  procedure  requires  sustained  load  testing  of  a  series  of  compact  specimens 
which  have  been  previously  fatigue  precracked.  Test  loads  for  the  individual  specimens  are 
selected  to  produce  a  range  of  values  of  initial  applied  stress  intensity  for  the  collection  of 
specimens.  Assuming  that  cracking  of  all  of  the  specimens  follows  a  mutual  propagation 
(da/dt,  K)  curve,  data  from  the  collection  of  tests  may  be  used  to  define  this  propagation 
curve.  The  data  required  include:  times  of  specimen  failure,  test  load,  initial  (a,)  and  final  (af) 
crack  length,  and  specimen  geometry.  The  following  mathematical  procedure  is  employed  to 
produce  the  crack  growth  curve: 

•  Initialize  the  coefficients  of  the  crack  growth  curve  at  values  which  are 
characteristic  of  the  material  in  question. 

•  Integrate  the  reciprocal  of  the  crack  growth  curve  from  initial  to  final 
crack  length  to  obtain  calculated  lives  (time  to  failure)  for  each  of  the 
individual  test  specimens. 

•  Calculate  the  summed  squared  error  for  the  actual  vs  observed  lives  for 
the  test  specimens. 

•  Generate  a  new  set  of  SINH  coefficients  using  the  Newton-Raphson  to 
minimize  the  calculated  summed  squared  error. 

•  Repeat  the  last  three  steps  above  until  convergence  is  achieved. 

This  procedure  was  applied  to  determine  the  sustained  load  crack  growth  behavior  of 
Waspaloy  (PWA  1007).  Nine  compact  specimens  were  tested,  five  at  649°C  (1200°F)  and  four 
at  732°C  ,  to  generate  the  crack  growth  model  shown  in  Figure  38.  (See  also  Figure  39.)  The 
range  of  validity  of  these  curves  is  approximately  45  MPa  \/~m  <  K<  110  MPa  \/~m,  and 
plotted  crack  growth  rates  are  approximately  equivalent.  The  crossing  of  the  curves  is  an 
artifact  of  the  range  of  testing,  and  extrapolation  beyond  the  range  of  the  data  should  not  be 
attempted. 

A  plot  of  the  actual  vs  calculated  times  to  failure  for  all  specimens  used  to  generate  the 
SINH  curves  is  presented  in  Figure  40.  Excellent  correlation  between  actual  and  calculated 
lives  is  shown  for  the  649°C  data.  However,  considerable  difficulty  was  encountered  in 
modeling  sustained  load  crack  propagation  at  732°C.  The  data  at  this  temperature  display 
much  less  consistency  than  at  649°C.  The  basis  for  the  inconsistency  probably  lies  with  a 
combination  of  factors  which  include  ill-defined  crack  fronts  and  crack  tunneling. 

(4)  Synergistic  Crack  Propagation 

(a)  LCF  —  Overload  Interaction 

Variable  amplitude  load  sequences  produce  synergistic  effects  on  crack  propagation 
which  complicate  the  life  prediction  procedure.  The  nature  and  significance  of  crack  retarda¬ 
tion  experienced  in  turbine  disks  are  discussed  in  Appendix  A.  Individual  SINH  models  des¬ 
cribing  the  effects  on  crack  growth  under  representative  load  sequences  are  presented  below. 
The  analysis  was  performed  on  data  generated  using  repetitive  mission  cycles  as  illustrated 
in  Figure  14.  As  discussed  previously  (and  in  Appendix  A),  the  overload-LCF  sequence  is 
defined  as  an  “elemental  damage  event.” 
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Figure  38.  Waspaloy  ( PWA  1007)  Crack  Propagation  Under  Sustained 
loading 
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Waspaloy  (PWA  1007),  Actual  vs  Calculated 


(1)  Effect  of  Number  of  Cycles  Between  Overloads 

The  number  (ANol)  of  baseline  fatigue  cycles  between  successive  overloads  has  a  pro¬ 
nounced  effect  on  the  average  crack  growth  rate  produced  by  the  repetitive  loading  of  Figure 
14.  For  tests  under  representative  load  sequence,  da/dN  vs  AK  data  were  generated  with 
da/dN  defined  as  the  average  crack  growth  per  cycle  (.iamiMion/(ANOL  +1))  and  AK  defined  for 
the  load  range  (Pmm  —  Pmin).  Thus,  the  repetitive  overload  is  treated  as  an  isolated  variable 
influencing  crack  growth,  and  the  AK  associated  with  the  overhead  is,  by  convention,  equiv¬ 
alent  to  the  baseline  AK. 

The  results  of  testing  of  Waspaloy  (PWA  1007)  conducted  under  conditions  of  R  =  0.5, 
0.167  Hz,  649°C  (1200°F),  OLR  =  POL/Pm  =  150'  ant*  5  —  AN0L  —  40  are  presented  in  a 
composite,  Figure  41.  Figure  42  lists  the  statistics.  Under  the  stated  conditions,  crack  growth 
rate  decreases  as  AN0L  increases.  The  SINH  curves  represent  an  interpolative  model  of  the 
effect  of  the  number  of  cycles  between  overloads.  As  shown  in  Figure  43,  each  of  the  defining 
SINH  coefficients  is  a  linear  function  of  (AN0L  +1). 

(2)  Effect  of  Overload  Ratio 

Waspaloy  (PWA  1007)  crack  propagation  data  generated  under  periodic  overload-fatigue, 
Figure  14,  with  overload  ratios  of  1.25  and  1.50  were  investigated.  The  other  conditions  of 
testing  were  equivalent:  ANol  =  40,  R  =  0.50,  0.167  Hz,  and  649  °C  U200°F).  Under  these 
conditions,  average  crack  growth  rate  increases  as  OLR  approaches  one.  When  OLR  is  equal 
to  one,  Kol  equals  KmM,  and  an  observable  ANol  vanishes.  Therefore,  for  OLR=  1.0  the  effect 
of  overload  ratio  reduces  to  baseline  crack  growth.  This  is  shown  in  Figure  44,  and  the  statis¬ 
tics  are  summarized  in  Figure  45.  The  SINH  curve  for  OLR  =  1.0  is  that  predicted  by  the 
constant  load  amplitude  stress  ratio  model  presented  earlier.  As  shown  in  Figure  46,  the 
SINH  coefficients  describing  the  effect  of  overload  ratio  are  a  linear  function  of  OLR. 

(b)  LCF  —  Dwell  Interaction 

Investigation  of  possible  synergistic  interaction  within  an  LCF-dwell  sequence  was  per¬ 
formed  using  the  repetitive  mission  shown  in  Figure  13.  Crack  propagation  tests  employing 
this  mission  were  conducted  for  Waspaloy  (PWA  1007)  at  649°C  (1 200°F)  and  732°C  (1350°F). 
The  crack  growth  data  were  reduced  with  da/dN  calculated  as  the  average  crack  advance 
per  cycle,  with  the  dwell  cycle  counted  as  one  load  cycle. 

The  results  of  three  649°C  (1200°F)  tests  of  10,  20,  and  40  sawtooth  fatigue  cycles 
between  dwells  of  120  sec  produced  similar  crack  growth  rates,  Figure  47.  Furthermore,  all  of 
this  data  is  adequately  described  by  the  previously  developed  model  of  baseline  fatigue  (0.167 
Hz,  R  =  0.10, 649°C).  Apparently,  there  is  little  interaction  or  effect  on  crack  growth  produced 
by  the  periodic  application  of  the  load  dwell  at  maximum  load  under  the  conditions  of 
testing. 

The  FCP  rate  under  similar  loading  at  732°C  (1350°F)  was  also  found  to  be  independent 
of  the  number  of  sawtooth  cycles  separating  the  periodic  load  dwells  of  120  sec.  The  data 
presented  in  Figure  48  represent  three  different  test  conditions:  10,  20,  and  40  cycles  between 
dwells.  All  of  this  data  is  described  by  the  single  SINH  curve  shown.  An  interpolative  model 
of  the  effect  of  temperature  on  crack  growth  under  LCF-Dwell  loading  is  illustrated  in  Figure 
49  and  the  statistics  shown  on  Figure  50.  As  shown,  crack  growth  rate  is  independent  of  the 
number  of  cycles  between  dwells  (10  <  AN^^,  <  40)  and  is  a  linear  function  of  temperature 
as  shown  in  Figure  51. 
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Figure  41.  Waspaloy  (PWA  1007)  Crack  Propagation,  Model  of  Effect  of 
the  Number  of  Cycles  Between  Overloads,  R  =  0.50,  0.167  Hz 
(10  cpm ),  649° C  (1200° F) 
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b.  GATORIZED-  IN100  (PWA  1073)  Crack  Propasatlon 

The  majority  of  the  IN100  crack  growth  data  presented  in  this  section  were  generated 
under  a  previous  AFML  contract  (References  1,  2,  and  3).  Additional  data  have  been  pro¬ 
duced  in  testing  under  the  current  contract,  and  data  from  both  sources  have  been  employed 
in  model  development.  The  models  are  presented  in  the  following  sections. 

(1)  Basic  Propagation,  Constant  Load  Amplitude  Cycling 

(a)  Effect  of  Loading  Rate  (Frequency) 

The  influence  of  cyclic  frequency  on  crack  growth  in  IN100  was  examined  over  the  range 
0.00833  (0.5  cpm)  <  frequency  <  20  Hz,  at  649°C  and  R  =  0.1.  An  interpolative  model  of  the 
effect  of  frequency  developed  from  this  data  is  shown  in  Figure  52  and  Figure  53  shows  the 
statistics.  The  SINH  coefficients  C2  and  C4  are  linear  functions  of  log  (frequency)  over  the 
full  range  of  the  model.  However,  this  is  not  the  case  for  the  coefficient  C3  which  locates  the 
abscissa  of  the  SINH  inflection  points.  As  illustrated  in  Figure  54,  C3  is  represented  as  a 
linear  function  for  0.00833  Hz  (0.5  cpm)  <  frequency  <  0.167  Hz  (10  cpm),  and  a  second 
function  is  defined  over  the  range  0.167  Hz  (10  cpm)  <  frequency  <  20  Hz.  This  relationship 
was  required  in  order  to  produce  an  adequate  description  of  the  data  over  the  range  of  test¬ 
ing.  The  0.00833  Hz  (0.5  cpm)  data  exhibit  anomalous  test-to-test  variability.  All  specimens 
were  tested  at  identical  conditions  (0.00833  Hz  (10  cpm),  R  =  0.10,  649°C);  however,  specimen 
604  AF  exhibited  a  much  slower  crack  growth  rate  than  the  population  of  tests.  In  fact,  the 
data  of  604  AF  agree  well  with  the  results  of  0.167  (10  cpm).  This  variability  cannot  be 
attributed  to  heat-to-heat  effects,  since  specimen  605  AF  (which  is  in  agreement  with  the 
remaining  data)  was  taken  from  the  same  heat  of  material  as  604  AF.  While  no  error  in  the 
604  AF  data  analysis  or  the  recorded  test  conditions  can  be  found,  an  unknown  error  in 
testing  is  believed  to  be  responsible  for  the  anomalous  behavior. 

(b)  Effect  of  Stress  Ratio,  0.167  Hz 

An  interpolative  model  of  the  effect  of  stress  ratio  on  crack  propagation  in  IN100  (PWA 
1073)  was  developed  from  data  over  the  range  0.1  <RS  0.8.  The  cyclic  frequency  was  0.167 
Hz  and  the  test  temperature  was  649°C.  The  base  condition  of  R  =  0.10  is  defined  by  a  collec¬ 
tion  of  data  from  six  tests  and  three  heats  of  material.  The  composite  SINH  model  for  the 
effect  of  stress  ratio  is  illustrated  in  Figure  55  and  Figure  56  gives  the  statistics.  The  asso¬ 
ciated  correlative  parameters  are  linear  functions  of  log  (1-R)  for  positive  stress  ratios  as 
shown  in  Figure  57. 

The  effect  of  negative  stress  ratio  fatigue  on  crack  propagation  of  IN100  is  discussed 
under  auxiliary  investigations. 

(c)  Effect  of  Stress  Ratio,  20  Hz 

The  effect  of  stress  ratio  on  crack  growth  in  IN100  (PWA  1073)  fatigued  at  a  frequency  of 
20  Hz,  is  described  by  the  model  shown  in  Figure  58,  and  Figure  59  reviews  the  statistics. 
The  defining  SINH  coefficients  are  illustrated  in  Figure  60. 

(d)  Effect  of  Temperature 

Results  of  crack  propagation  testing  of  IN100  (PWA  1073)  conducted  at  0.167  Hz,  R  =  0.1, 
and  538°C  (1000°F)  STS  732°C  (1350°F),  were  used  to  develop  an  interpolative  SINH  model 
of  the  effect  of  temperature.  This  model  is  presented  in  Figure  61,  the  statistics  presented  in 
Figure  62,  and  the  defining  SINH  coefficients  are  given  in  Figure  63. 
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Figure  53.  Statistics  for  IN  100  ( PWA  1073)  Crack  Propagation,  Frequency  Model,  R  =  0.10,  649° C  (1200°F) 


FO  169787 

Figure  54.  IN  100  (PWA  1073)  Crack  Propagation,  Frequency  Model  Correl- 
j  ative  Parameters,  R  =  0.10,  649° C  (1200°F) 
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Figure  59.  Statistics  for  IN  100  ( PWA  1073)  Crack  Propagation.  Stress  Patio  Model.  20  Hz.  OV>°C  (1 200" l  ) 
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Figure  62.  Statistics  for  IN 100  (PW A  1073)  Crack  Propagation,  Temperature  Model,  U  -  0.10,  0.167  llz  (10  cpm) 


COEFFICIENT 


(2)  Crack  Propagation  Under  Periodic  Load  Dwell 

(a)  Effect  of  Load  Dwell,  649°C  (1200°F) 

Crack  propagation  under  periodic  load  dwell  in  GATORIZED®  INlOO  (PWA  1073)  was 
investigated  using  the  repetitive  loading  shown  in  Figure  12.  The  period  of  load  dwell  was 
60, 120,  and  600  sec,  in  three  tests  where  R  =  0.10  and  T  =  649°C.  The  composite  model  of  the 
effect  of  the  period  of  load  dwell  is  illustrated  in  Figure  65.  Figure  64  lists  the  statistics.  The 
interpolative  SINH  coefficients  are  shown  in  Figure  66. 

(b)  Effect  of  Load  Dwell,  732°C  (1350°F) 

The  effect  of  period  of  load  dwell  on  crack  growth  in  INlOO  tested  at  732°C  is  described 
by  the  interpolative  SINH  model  presented  in  Figure  67.  Figure  68  gives  the  statistics.  The 
increased  crack  growth  rate  associated  with  increased  dwell  period  is  correlated  with  the 
SINH  coefficients  as  shown  in  Figure  69. 

(3)  Sustained  Load  Crack  Growth 

Under  a  previous  Air  Force  program  the  sustained  load  crack  growth  behavior  of  INlOO 
(PWA  1073)  was  determined  (Reference  2).  The  experimental  and  analytical  procedures  used 
were  the  same  as  described  earlier  in  this  report  for  Waspaloy  sustained  load  crack  growth. 
The  results  of  the  INlOO  investigation,  presented  in  Figures  70  and  71,  demonstrate  the 
deleterious  influence  of  increasing  temperature  on  sustained  load  crack  growth.  The  correla¬ 
tion  between  these  curves  and  10  tests  of  12.7  mm  (0.50  in.)  thick  compact  specimens  is 
illustrated  in  Figure  72. 

(4)  Synergistic  Crack  Propagation 
(a)  LCF  —  Overload  Interaction 

The  nature  of  crack  retardation  resulting  from  load  sequencing  representative  of  mil¬ 
itary  gas  turbine  operation  is  discussed  in  Appendix  A.  Testing  based  on  the  repetitive  load 
cycle  shown  in  Figure  14  was  used  to  characterize  synergistic  interaction  for  LCF-overload 
fatigue  of  GATORIZED®  INlOO.  Interpolative  models  of  crack  propagation  as  a  function  of 
overload  occurence  frequency  and  relative  magnitude  are  presented  below. 

(1)  Effect  of  Cycles  Between  Overloads 

Under  constant  load  amplitude  fatigue  (R  =  0.50)  at  649°C  (1200°F)  periodically  inter¬ 
rupted  by  a  single  overload  (OLR  =  1.5),  crack  propagation  in  INlOO  displays  a  strong 
dependence  on  the  number  of  baseline  fatigue  cycles  between  overloads  (AN0L).  Results  of 
crack  propagation  tests  conducted  with  ANol  =  5,  20,  and  40  cycles  are  plotted  in  composite 
fashion  in  Figure  73.  Figure  74  presents  the  statistics.  The  associated  SINH  curves  represent 
an  interpolative  model  of  the  effect  of  number  of  cycles  between  successive  periodic  over¬ 
loads.  As  shown  in  Figure  75,  the  defining  SINH  coefficients  are  linear  functions  of  the 
logarithm  of  the  number  of  cycles  per  repetitive  mission  (log(ANOL  +  1)). 
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IN100  (PWA  1073)  Crack  Propagation,  Correlative  Parameters 
Model  of  Effect  of  £sNQL,  R  =  0.50,  0.167  Hz  (10  cpm),  649° C 
(1200° F).  OLR  =  1.50 
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(2)  Effect  of  Overload  Ratio 

Under  periodic  overload-fatigue  the  effect  of  increasing  overload  ratio  (POL/PraM)  while 
holding  ANq[  constant  is  to  decrease  the  average  rate  of  crack  propagation  produced  by  the 
load  sequencing.  Tests  of  IN  100  were  conducted  using  the  repetitive  load  sequence  shown  in 
Figure  14.  Similar  to  the  tests  discussed  above,  the  temperature  was  649°C  (1200°F),  R  =  0.50, 
and  the  cyclic  frequency  was  0.167  Hz  (10  cpm).  The  number  of  cycles  between  overloads, 
AN0L,  was  held  constant  at  40,  and  overload  ratio  was  varied  over  the  range  1.0  <  OLR  < 
1.50. 

Under  an  overload  ratio  of  1.0,  the  load  range  is  constant  and  crack  growth  reduces  to 
the  case  of  the  constant  AP  baseline.  For  overload  ratios  of  1.25  and  1.50  average  growth  rate 
is  reduced.  A  composite  of  these  data,  and  the  associated  SINH  model,  is  presented  in  Figure 
76.  Figure  77  gives  the  statistics.  The  defining  SINH  coefficients  are  linear  functions  of  OLR 
as  illustrated  in  Figure  78. 

( b )  LCF  —  Dwell  Interaction 

Crack  propagation  testing  of  IN100  fatigued  at  649°C  (1200°F)  under  the  repetitive  LCF- 
dwell  sequence  shown  in  Figure  13  revealed  little  or  no  effect  of  the  periodic  dwell.  As  illus¬ 
trated  in  Figure  79,  data  generated  with  10,  20,  and  40  cycles  of  0.167  Hz  (10  cpm)  sawtooth 
fatigue  (R  =  0.1)between  dwells  correlated  well  with  the  crack  growth  curve  representing 
sawtooth  fatigue  only.  Figure  80  presents  the  statistics. 

Similar  tests  conducted  at  732°C  (1350  °F)  demonstrated  a  mild  effect  of  the  periodic 
load  dwell.  The  results  of  tests  conducted  with  10,  20,  and  40  cycles  between  dwells  are 
shown  in  a  composite  presented  in  Figure  81.  Figure  82  presents  the  statistics.  For  purposes 
of  comparison,  the  SINH  curves  describing  crack  growth  under  baseline  fatigue  (732°C, 
0.167,  R  =  0.10)  and  repetitive  dwell  (732°C,  120  sec  dwell,  R  =  0.10)  are  also  shown.  The  SINH 
curves  representing  the  individual  LCF-dwell  data  sets  are  given  by  an  interpolative  model 
of  this  effect,  and  the  defining  coefficients  are  illustrated  in  Figure  83. 

5.  Auxiliary  Investigation 
a.  Negative  Strata  Ratio  Ettacta 

Tension-compression  fatigue  (R  <  O)  may  occur  in  isolated  disk  locations  during  mission 
operation  as  shown  in  Figure  4.  Fatigue  crack  propagation  under  such  negative  stress  ratio 
cycling  is  not  accurately  predicted  by  a  simple  extrapolation  of  a  stress  ratio  model  devel¬ 
oped  for  tension-tension  (R  >  O)  fatigue.  Due  to  a  difference  in  the  mechanism  of  crack 
growth  associated  with  the  positive  vs  negative  load  excursion,  a  separate  evaluation  of 
crack  growth  is  required  for  R  <  0. 

Considering  the  work  of  Elber  (Reference  7),  the  effective  stress  intensity  range  operat¬ 
ing  at  the  crack  tip  is  AK^  =  (Kmtx  •  Kclo,ore),  where  Kclo,ore  is  the  stress  intensity  at  which 
physical  closing  of  the  crack  tip  occurs  on  unloading.  Under  a  compressive  stress  field,  a 
crack  tip  singularity  does  not  exist,  and  K  =  0  for  all  <7,pplied  <  O.  While  for  a  given  Kmtt,  the 
applied  stress  intensity  range  (AK  =  Krau  -  K^)  remains  constant  for  negative  stress  ratio 
fatigue,  is  increased  slightly  as  decreases  due  to  compressive  yielding. 

In  order  to  examine  the  effect  of  representative  tension-compression  fatigue  on  crack 
growth  in  Waspaloy  and  IN100,  the  matrix  of  tests  outlined  in  Table  5  was  conducted.  The 
investigation  was  performed  at  427°C  (800°F)  and  649°C  (1200°F)  and  spanned  the  range  of 
-1.0  <  R  S  0.0;  the  cyclic  frequency  was  0.167  Hz  (10  cpm).  The  crack  growth  specimen 
employed  in  these  tests  was  of  through-thickness  center  crack  geometry,  Figure  11a. 
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Figure  H2.  Statistics  for  IN  100  (PWA  1073)  Crack  Propagation,  Model  of  LCF-Dwell  Interaction,  H  ~  0.10,  732° C  (1330° F) 
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All  tests  were  conducted  with  the  specimen  carefully  aligned  and  ends  rigidly  fixed.  In 
spite  of  such  efforts  to  achieve  a  uniform  nominal  stress  field  in  the  specimen,  crack  growth 
irregularities  observed  in  several  specimens  suggested  the  presence  of  a  bending  stress.  The 
irregular  crack  growth  was  characterized  by  unequal  rates  of  propagation  at  the  opposing 
specimen  surfaces,  producing  crack  front  asymmetry.  This  condition  required  the  rejection  of 
data  from  several  completed  tests  and  subsequent  retesting. 

(1)  Waspaloy 

The  results  of  the  427°C  (800°F)  negative  stress  ratio  tests  for  Waspaloy  are  presented  in 
Figure  84.  The  rate  of  crack  propagation  is  observed  to  be  similar  for  R  =  -0.5  and  R  =  -1.0 
fatigue,  and  crack  growth  under  both  of  these  stress  ratios  is  more  rapid  than  the  illustrated 
curve  for  R  =  0  cycling.  Since  all  data  are  regressed  with  Kmin  =  0  for  compressive  stress 
excursions,  the  comparison  of  the  negative  R  data  with  the  R  =  0  curve  isolates  the  influence 
of  the  compressive  excursion. 

The  effect  of  negative  stress  ratio  cycling  on  crack  propagation  of  Waspaloy  at  649°C 
(1200°F)  is  illustrated  in  Figure  85.  Again  the  compressive  loading  is  observed  to  increase 
crack  growth  above  the  rate  produced  by  R  =  0  fatigue. 

(2)  IN  100 

Crack  propagation  in  IN100  subject  to  negative  stress  ratio  exhibited  anomalous  behav¬ 
ior.  At  427°C  (800°F),  Figure  86,  the  effect  of  R  =  -0.5  fatigue  is  nearly  equivalent  to  behavior 
under  R  =  0  cycling,  however,  decreasing  stress  ratio  to  R  =  -1.0  (specimen  793)  is  observed  to 
reduce  crack  growth  rate  significantly  below  that  of  baseline  (R  =  0)  fatigue.  Close  examina¬ 
tion  of  the  data  from  specimen  793  revealed  no  satisfactory  explanation  for  the  slower  crack 
growth.  While  this  specimen  was  taken  from  a  separate  heat  of  IN  100,  experience  with  this 
alloy  indicates  that  the  anomalous  behavior  cannot  be  explained  as  a  heat-to-heat  effect. 

Anomalous  crack  growth  was  also  observed  in  649°C  (1200°F)  negative  stress  ratio  test¬ 
ing  of  IN100.  As  illustrated  in  Figure  87,  crack  propagation  under  R  =  -0.5  fatigue  was 
observed  to  be  more  rapid  than  under  R  =  -1.0  cycling.  Propagation  under  both  conditions 
was  more  severe  than  for  R  =  0  fatigue. 

(3)  Interpol  alive  Modeling  of  Crack  Growth 

The  basis  for  the  anomalies  in  crack  propagation  exhibited  by  the  negative  stress  ratio 
data  is  believed  to  be  the  tendency  for  asymmetric  crack  growth,  as  previously  discussed. 
While  the  data  presented  were  taken  from  specimens  displaying  very  limited  crack  asymme¬ 
try,  such  behavior  appears  to  be  the  source  of  appreciable  test-to-test  variability  for  tension- 
compression  fatigue.  Similar  crack  asymmetry  and  data  variability  were  not  observed  under 
tension-tension  fatigue.  The  use  of  a  large  capacity  testing  machine,  having  a  very  rigid  load 
train,  should  significantly  reduce  the  variability  for  tension-compression  fatigue.  Such 
equipment  was  not  available  for  the  current  program. 

In  the  light  of  the  observed  variability  in  the  generated  crack  propagation  data  for  nega¬ 
tive  stress  ratio  fatigue,  interpolative  models  of  crack  growth  behavior  could  not  be  reliably 
developed.  Multiple  tests  at  each  condition  are  suggested  in  order  to  appropriately  character¬ 
ize  the  effect  of  tension-compression  cycling. 
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Figure  84.  Waspaloy  (PWA  1007)  Crack  Propagation,  Negative  Stress  Ratio 
Effect,  0.167  Hz  (10  cpm),  427° C  (800° F) 
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b.  E fleet  of  Net  Section  Strete 


The  effect  of  high  net  section  stress  on  crack  propagation  in  Waspaloy  (PWA  1007)  was 
investigated  using  the  through-thickness  center  crack  specimen  in  Figure  1  la.  This  specimen 
was  selected  for  its  small  obtainable  crack  size  and  the  ability  to  achieve  high  net  section 
stresses  at  stress  intensities  comparable  with  those  of  the  baseline  testing.  Tests  of  two  spec¬ 
imens  were  conducted  with  onn  >  0.8  oylrtd,  T  =  649°C  (1200°F),  R  =  0.05,  and  a  frequency  of 
0.167  Hz  (10  cpm).  The  resulting  data  are  presented  in  Figure  88,  and  their  relationship  to  the 
baseline  crack  growth  curve  is  illustrated.  Crack  propagation  under  the  high  stress  loading 
is  observed  to  be  slightly  less  rapid  than  under  equivalent  test  conditions  with  significantly 
reduced  stress.  Thus,  the  interpolative  models  of  crack  propagation  in  Waspaloy  (PWA  1007) 
should  be  applicable,  but  slightly  conservative  when  used  to  predict  crack  growth  in  a  stress 
field  of  magnitude  approaching  the  material  yield  strength. 

Similar  findings  have  been  reported  (Reference  2)  for  the  effect  of  high  net  section  stress 
on  elevated  temperature  crack  propagation  in  GATORIZED?  IN100  (PWA  1073). 

c.  Prior  Plettlc  Detormetlon 

Component  crack  propagation  life  predictions  are  made  using  da/dN  data  obtained  from 
virgin  (unstrained)  material.  However,  cracks  occuring  in  real  hardware  often  initiate  in 
regions  of  high  stress  concentration  such  as  boltholes  or  blade  attachment  areas.  Material  in 
these  regions  may  be  subjected  to  many  strain  cycles  in  excess  of  1.0%.  It  is  assumed  that 
FCP  behavior  (da/DN  vs  aK)  in  this  material  is  essentially  the  same  as  that  for  virgin 
material  and  that  data  generated  on  virgin  specimens  can  be  used  to  describe  FCP  in  pre¬ 
strained  material.  The  purpose  of  this  auxiliary  investigation  was  to  test  this  hypothesis, 
using  Waspaloy  (PWA  1007)  and  IN100  (PWA  1073). 

(1)  Waspaloy 

In  order  to  simulate  the  cyclic  loading  experienced  in  the  vicinity  of  a  bolthole  in  a 
Waspaloy  turbine  disk,  1.0%  tensile  strain  range  cycling  (t  meM  =  %  tmu)  was  applied  to 
unnotched  center-crack  specimens  (Figure  11a).  The  specimens  (1321  and  1323)  were  cycled 
at  649°C  (1200°F)  and  0.167  Hz  (10  cpm)  for  10%  (200  cycles)  of  the  expected  crack  initiation 
life.  Determination  of  the  cyclic  stress-strain  behavior  of  Waspaloy  was  performed  using  the 
strain  control  specimen  of  Figure  89,  and  typical  hysteresis  loops  tire  presented  in  Figure  90. 
During  the  elevated  temperature  pre-strain  cycling,  specimen  1321  was  monitored  with 
acoustic  emission  equipment,  and  a  plot  of  the  resulting  emission  count  is  presented  in  Fig¬ 
ure  91.  The  initially  high  acoustic  activity  attenuated  to  a  stable  level  within  the  first  50 
cycles.  Acoustic  emission  reflects  damage  accumulation  during  testing.  Ho  detectable  flaws 
resulted  from  the  pre-straining. 

The  unnotched  specimens  were  removed  from  the  test  machine,  and  a  through-thickness 
center  flaw  was  machined  in  each  using  electric  discharge  methods.  From  this  point  the 
specimens  were  treated  the  same  as  virgin  specimens.  The  notched  specimens  were  returned 
to  the  testing  machine,  precracked  at  room  temperature,  and  the  crack  was  propagated  to 
failure  under  constant  load  amplitude  fatigue  T  =  649°C  (1200°F),  0.167  Hz  (10  cpm),  and  a 
stress  ratio  of  0.1.  The  crack  propagation  data  (da/dN  vs  AK)  and  a  SINH  curve  representing 
crack  propagation  in  virign  (unstrained)  Waspaloy  are  presented  in  Figure  92.  Crack  growth 
in  the  pre-strained  material  is  approximately  30%  faster  than  in  virgin  material.  The  crack 
growth  rate  in  the  pre-strained  specimens  is  approximately  30%  faster  than  observed  from 
the  virgin  material  tests.  However,  this  deviation  is  explainable  by  intrinsic  data  scatter  for 
crack  propagation  testing. 
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Figure  88.  Waspaloy  (PWA  1007)  Crack  Propagation,  High  Stress  Effect, 

R  =  0.05,  0.167  Hz  (10  cpm),  649° C  (1200° F) 
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Figure  89.  Strain  Controlled  Low  Cycle  Fatigue  Specimen 


Figure  90.  Hysteresis  Observed  in  1.0%  Strain  Controlled  Testing 


Acoustic  Emission  Rate  (Counts/Cycle) 


Figure  91.  Acoustic  Emission  Record  of  Prestrain  Cycling  of  Waspaloy 
(PWA  1007)  Specimen 
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Figure  92.  Waspaloy  (PWA  1007)  Crack  Propagation,  Effect  of  Prior . 

Plastic  Deformation,  R  =  0.10,  0.167  Hz  (10  c pm),  649° C 
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(2)  IN100 


The  effect  of  1.0%  strain  cycling  prior  to  crack  propagation  was  also  examined  in  IN100. 
The  procedure  was  as  described  for  Waspaloy,  and  the  results  are  illustrated  in  Figure  93. 
Crack  propagation  in  pre-strained  IN  100  was  observed  to  be  generally  equivalent  to  charac¬ 
teristic  crack  growth  in  virgin  IN100. 

(3)  Component  Lite  Prediction 

The  results  of  the  Waspaloy  and  IN100  tests  discussed  above  indicate  different  effects  of 
1.0%  pre-strain  cycling.  Using  a  crack  growth  model  developed  for  virgin  material  data,  pre¬ 
diction  of  fatigue  crack  life  in  pre-strained  Waspaloy  (PWA  1007)  are  expected  to  be  slightly 
anticonservative,  while  for  IN100,  little  model  bias  is  indicated.  However,  the  tests  results  for 
both  materials  fall  within  the  range  of  intrinsic  data  scatter.  Therefore,  fatigue  crack  growth 
models  representing  virgin  material  data  from  both  alloys  are  applicable  to  life  prediction  in 
pre-strained  material. 

d.  Overload-Underload-L  CF  Interaction 

Flight  mission  stress  analysis  which  was  conducted  under  Phase  I  of  this  contract 
revealed  the  existence  of  compressive  stresses  within  the  loading  spectrum  of  a  cooled  tur¬ 
bine  disk  (Figure  4).  The  compressive  stresses  occurred  immediately  following  a  major  tensile 
load  excursion  (overload)  and  were  expected  to  have  a  detrimental  influence  upon  crack 
propagation  rate  by  reducing  the  beneficial  effect  of  the  overload. 

A  test  was  conducted  to  examine  the  effect  on  crack  propagation  of  compressive  loading 
in  the  repetitive  mission  sequence  of  overload-underload-low  cycle  fatigue,  Figure  94.  The 
resulting  data  was  reduced  with  the  hyperbolic  sine  model  and  is  presented  in  the  compara¬ 
tive  plot  of  Figure  95.  The  three  data  sets  are  from  tests  of  IN100  at  649°C  (I200°F)  and 
(1)  constant  load  amplitude  fatigue,  R  =  0.5;  (2)  constant  load  amplitude  fatigue,  R  =  0.5, 
interrupted  by  a  periodic  overload  (Overload  Ratio  =  P0[/Pm„  =  1.5)  every  21  cycles;  and 
(3)  identical  to  (2)  with  the  addition  of  a  single  compressive  load  excursion  with  an  underload 

of  (Pui/P*  J=‘ -0-45- 

Crack  propagation  resulting  from  the  overload  LCF  loading  (curve  2)  is  significantly 
less  rapid  than  under  similar  constant  load  amplitude  fatigue  (curve  1).  However,  the 
overload-underload-LCF  mission  produces  crack  propagation  which  is  slightly  more  severe 
than  under  constant  AP  fatigue.  That  is,  in  IN  100  (PWA  1073)  cycled  at  649°C  (1200°F)  the 
effect  of  the  minor  compressive  load  excursion  (U0L  =  -0.45)  following  an  overload  (OLR  =  1.5, 
ANol  =  20)  is  to  erase  the  beneficial  effect  of  the  overload,  and  the  combined  effect  of  the 
overload-underload  sequence  is  detrimental  to  the  crack  propagation  life. 

Further  examination  of  the  influence  of  overload-underload-LCF  sequencing  on  crack 
growth  is  beyond  the  scope  of  the  present  contract.  Such  research  is  suggested  for  a  future 
program. 

a.  Thermal-Mechanical  Fatigue 

Because  the  study  of  thermal-mechanical  synergism  is  particularly  complex,  it  has  been 
assumed  that  the  effects  of  frequency,  stress  ratio,  dwells,  and  overloads  occurring  with 
changing  temperature  can  be  approximated  by  considering  the  isothermal  increments  of 
each  parameter.  The  results  of  a  survey  test,  discussed  below,  tend  to  substantiate  this 
approach. 
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Figure  93.  IN  100  (PWA  1073)  Crack  Propagation,  Effect  of  Prior  Plastic 
Deformation,  R  =  0.10,  0.167  Hz  (10  cpm),  649° C  (1200° F) 
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The  thermal-mechanical  fatigue  (TMF)  test  conducted  under  this  task  employed  saw¬ 
tooth  waveforms  for  both  thermal  and  mechanical  cycling.  Temperature  and  load  cycling 
were  applied  180°  out-of-phase  as  shown  in  Figure  96,  and  the  test  specimen  was  of  through¬ 
thickness  center  crack  geometry  as  shown  in  Figure  11a.  The  applied  cyclic  frequency  was 
0.0167  Hz  (1  cpm)  which  required  heating  and  cooling  rates  of  556°C  (1000CF)  per  minute  for 
temperature  cycling  between  427°C  (800°F)  and  704  °C  (1300°F). 
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Figure  96.  Thermal-Mechanical  Fatigue  Cycle 


Specimen  heating  was  accomplished  with  electrical  induction,  and  forced  air  was  used  to 
achieve  the  proper  cooling  rate.  An  investigation  of  the  through-thickness  thermal  gradient 
in  the  specimen  was  conducted  by  mounting  thermocouples  on  the  specimen  surface  and  in  a 
drill  hole  located  at  the  specimen  center.  The  internal  temperature  lagged  the  surface 
temperature  by  a  maximum  of  56°C  (100°F)  for  a  556°C  (1000°F)  per  minute  cooling  rate. 

Figure  97  compares  the  thermal-mechanical  FCP  observed  under  out-of-phase  cycling 
(maximum  temperature  occurring  at  minimum  load)  with  crack  growth  under  isothermal 
conditions  at  the  same  frequency  and  stress  ratio  (0.0167  Hz  (10  cpm),  R  =  0.10).  The  four  SINH 
curves  represent  isothermal  crack  growth  at  a  series  of  temperatures  which  spans  the  range  of 
thermal  cycling  (427°C  (800°F)  to  704  °C  (1300°F)).  The  TMF  data  are  observed  to  correlate  well 
with  results  of  isothermal  tests  conducted  at  427°C.  This  temperature  coincides  with  the 
temperature  at  which  peak  loading  occurs  during  the  TMF  cycle,  signifying  the  importance  of 
the  peak  load  vs  temperature  condition.  This  result  indicates  that  TMF  behavior  is  more  closely 
related  to  isothermal  FCP  than  originally  thought. 


1.  Etfact  of  Spmciman  Thlcknott 

An  investigation  of  the  effect  of  specimen  thickness  on  elevated  temperature  CTack  propa¬ 
gation  of  Waspaloy  (PWA  1007)  was  conducted  at  427°C  (800°F)  and  649°C  (1200°F).  At  each 
temperature,  specimens  of  three  thicknesses  (2.54,  7.62,  and  12.70  mm)  (0.1,  0.3,  and  0.5  in.) 
were  tested  at  similar  conditions  (0.0167  Hz  (10  cpm),  R  =  0.05).  The  results  of  these  tests  are 
illustrated  in  Figure  98  and  99.  At  each  temperature,  equivalent  rates  of  crack  growth  are 
observed  for  the  7.62  and  12.70  mm  thick  specimens,  while  a  much  slower  rate  of  crack 
propagation  occurred  in  the  2.54  mm  thick  specimen.  From  these  findings  it  is  concluded 
that  the  rate  of  FCP  of  Waspaloy  (PWA  1007)  becomes  thickness  independent  at  some  thick¬ 
ness  less  than  7.62  mm  (0.30  in.). 
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Figure  98.  Waspaloy  (PWA  1007)  Crack  Propagation,  Effect  of  Specimen 
Thickness,  R  =  0.05,  0.167 Hi  (10  cpm),  427° C  (800° F) 
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Figure  99.  Waspaloy  (PWA  1007)  Crack  Propagation,  Effect  of  Specimen 
Thickness,  R  =  0.05,  0.167  Hz  (10  cpm),  649° C  (1200°F) 
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D.  PHASE  IV  -  MODEL  DEMONSTRATION 


The  goal  of  this  phase  of  the  program  was  to  demonstrate  the  accuracy  of  the  developed 
crack  propagation  model.  During  the  initial  step  of  this  effort,  the  Air  Force  project  engineer 
provided  two  mission  profiles  upon  which  life  history  (a  vs  N)  predictions  were  made.  Subse¬ 
quently,  crack  propagation  tests  of  Waspaloy  (PWA  1007)  and  IN100  (PWA  1073)  were  con¬ 
ducted  under  repetitive  application  of  the  provided  missions  in  two  specimen  geometries  — 
compact  and  surface  flaw. 

1.  Demonstration  Missions 

The  two  demonstration  missions  supplied  by  the  Air  Force  are  illustrated  in  Figures  100 
and  101.  In  order  to  fully  exercise  the  predictive  capability  ot  the  crack  growth  model,  these 
missions  were  chosen  to  be  extreme  cases.  They  were  not  defined  to  simulate  enigne  opera¬ 
tion.  The  first  mission,  Figure  100,  was  composed  of  a  series  of  major  load  excursions,  fol¬ 
lowed  by  series  of  load  dwells  of  increasing  magnitude,  and  ending  with  a  flurry  of  fatigue 
cycles  at  a  high  load  ratio.  This  was  contrasted  with  the  second  demonstration  mission, 
shown  in  Figure  101,  which  contained  much  more  prominent  load  sequence  effects. 

2.  Mission  Segregation  and  Llfs  Prediction 

Cumulative  damage  segmentation  of  the  model  demonstration  missions  presented  in 
Figure  100  and  101  was  accomplished  by  use  of  the  segregation  algorithm  discussed  earlier. 
This  computerized  procedure  characterized  the  individual  loading  cycles  in  ternas  of  basic 
fatigue  parameters  such  as  maximum  load  (PmM),  load  ration  (r),  and  frequency.  Further¬ 
more,  the  segregation  algorithm  accounted  for  load  sequence  effects  which  result  in  crack 
growth  synergism.  Major  load  excursions  were  identified  and  load  sequence  parameters  des¬ 
cribing  overload  ratio  and  occurence  frequency  were  defined  by  the  algorithm. 

For  example,  model  demonstration  mission  number  1  begins  with  four  major  load  excur¬ 
sions.  These  were  defined  in  terms  of  Pm„,  R,  and  frequency,  and  the  fourth  of  these  cycles 
were  identified  as  an  overload  affecting  subsequent  cycles.  All  of  the  remaining  cycles  in  this 
mission  were  of  reduced  maximum  load  and  are  therefore  influenced  by  the  overload.  The 
number  of  overload  affected  cycles  (ANol)  was  defined,  and  the  overload  ratios 
were  calculated  for  the  individual  post-overload  cycles.  Note  that  the  subsequent  repetitive 
mission  marks  the  end  of  the  cycles  influenced  by  the  overload,  since  the  overload  ratio 
returns  to  1.0 

The  occurrence  of  the  series  of  stairstep  load  dwells  within  mission  1  posed  an  unusual 
situation.  As  discussed  earlier,  testing  under  LCF-Dwell  load  sequencing  revealed  that  the 
crack  growth  resulting  from  the  intermittent  load  dwells  was  significantly  slower  than 
observed  for  equivalent  repetitive  load  dwells.  This  behavior  occurred  in  both  Waspaloy  and 
IN100  and  agrees  with  the  findings  of  Macha,  Grandt,  and  Wicks.  (Reference  10.)  These 
researchers,  performing  elevated  temperature  crack  propagation  testing  of  IN100,  observed 
minimal  effect  of  load  dwells  in  a  representative  mission  load  sequence.  This  finding  is 
incorporated  in  the  segregation  algorithm  by  defining  as  zero  the  crack  growth  associated 
with  sustained  load  mission  segments.  Crack  propagation  produced  by  a  load  dwell  cycle 
was  assumed  to  result  entirely  from  the  initial  load  excursion  in  the  dwell  cycle.  In  the  event 
that  successive,  increasing  load  dwells  occurred,  the  stairstep  loading  sequence  was  reduced 
to  a  single  monotonic  load  excursion.  Therefore,  the  stairstep  sequence  of  mission  1  was 
characterized  as  a  single  loading  cycle  which  occurred  during  the  indicated  time  interval. 
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Figure  100.  Model  Demonstra:  i>n  Mission,  Xo. 


Segregation  of  model  demonstration  mission  number  2  (Figure  101)  was  also  performed 
by  the  computer  algorithm.  All  cycles  were  defined  in  terms  of  load,  load  ratio,  and  fre¬ 
quency,  and  the  first  cycle  was  identified  as  an  overload.  The  lsat  cycle  in  the  mission,  being 
equivalent  to  the  overload  cycle,  terminated  the  calculated  effect  of  the  initial  overload.  The 
number  of  overload  affected  crack  growth  cycles  was  determined  and  the  individual  ratios 
were  calculated. 

The  mission  segregations  discussed  above  were  the  basis  of  a  synergistic  model  life  pre¬ 
diction  accounting  for  load  sequence  interactions.  A  cycle-by-cycle  partitioning  of  the  test 
missions  assuming  that  load  sequence  does  not  affect  crack  growth  (i.e.  that  linear  damage 
summation  method)  was  performed  by  truncating  the  load  sequence  parameters  from  the 
output  of  the  synergistic  segregation  algorithm.  This  abbreviated  mission  segregation 
formed  the  basis  of  the  SINH  Linear  Damage  Model  life  prediction. 

The  segregated  missions  were  subsequently  input  into  the  interpolation  algorithm  which 
defined  a  characteristic  crack  growth  curve  for  each  unique  loading  condition.  The  output  of 
this  routine  was  that  used  in  a  cycle-by-cycle  integration  of  crack  growth,  known  as  the 
Computational  Algorithm,  to  obtain  life  calculations  for  a  total  of  eight  demonstration  speci¬ 
mens.  There  were  four  specimens  in  each  material,  Waspaloy  (PWA  1007)  and  IN100  (PWA 
1073).  Two  specimen  geometries  (compact  and  surface  flaw)  were  tested  under  each  of  the  two 
missions  (Figures  100  and  101).  Since  a  linear  damage  and  a  synergistic  model  life  prediction 
was  calculated  for  each  of  the  eight  specimens,  a  total  of  sixteen  life  calculations  were  per¬ 
formed.  These  were  recorded  with  the  Air  Force  project  engineer  prior  to  mechanical  testing. 


3.  Demonstration  Tasting 

The  two  specimen  geometries  used  in  demonstration  testing  were  the  compact  (Figure 
lib)  and  the  surface  flaw  specimen  (Figure  102).  The  compact  specimen  has  a  well  character¬ 
ized  stress  intensity  solution  (Reference  5)  and  was  used  extensively  during  the  test  program, 
while  the  surface  flaw  specimen  is  less  familiar  and  was  unused  in  data  generation.  The 
stress  intensity  solution  for  this  specimen  was  taken  from  calculations  performed  by  P&WA/ 
Commercial  Products  Division  (Reference  11).  Figure  103  compares  this  analytically  deter¬ 
mined  surface  flaw  K-solution  and  a  solution  calculated  using  handbook  values  (Reference 
12). 


Experimental  verification  of  the  accuracy  of  the  stress  intensity  solution  was  performed 
by  testing  a  Waspaloy  (PWA  1007)  surface  flaw  specimen  at  fatigue  conditions  equivalent  to 
baseline  data  generated,  under  the  test  program  of  this  contract.  This  data  is  illustrated  in 
Figure  104  with  2b/w  <  0.45,  and  good  correlation  with  the  SINH  curve  generated  from 
baseline  data  is  observed. 

All  testing  was  performed  on  servocontrolled  hydraulic  machines  operated  under  load 
control  with  mission  profile  described  by  either  an  MTS  minicomputer  or  a  DATATRAK® 
controller.  The  atmosphere  was  laboratory  air,  and  specimen  heating  was  accomplished 
using  clamshell  resistance  furnaces. 
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Figure  103.  Comparison  of  K-Calibration  Curves  for  Surface  Flaw  Specimen 


4.  Model  Verification 

Model  demonstration  tests  were  conducted  to  determine  the  accuracy  of  the  crack  propa¬ 
gation  models  developed  under  this  program.  The  procedure  of  modeling  crack  growth  under 
complex  loading  conditions  was  equivalent  for  both  Waspaloy  and  IN  100.  The  test  data  from 
all  specimens  were  treated  in  a  uniform  fashion.  In  order  to  eliminate  transient  effects  on 
crack  propagation  commonly  associated  with  change  from  room  temperature  precracking 
conditions  to  elevated  temperature  test  conditions,  the  initial  0.635  mm  (0.025  in.)  of  crack 
growth  was  deleted  from  all  data  sets. 

s.  Watpaloy  (PWA  1007) 

All  mission  testing  to  evaluate  the  effectiveness  of  the  Waspaloy  (PWA  1007)  crack 
propagation  model  was  conducted  at  a  single  temperature:  621°C  (1150°F).  The  results  of 
tests  performed  under  Mission  1  (Figure  100)  using  compact  and  surface  flaw  specimens  are 
presented  in  Figures  105  and  106,  respectively.  For  both  specimens,  the  life  predictions  agree 
well  with  the  data. 

A  parameter  which  may  be  used  to  sample  the  accuracy  of  the  live  prediction  is  the  ratio 
of  calculated  cycles  to  failure  vs  actual  cycles  to  failure  (NC/NA).  For  the  compact  specimen 
(Figure  105)  this  parameter  demonstrates  the  superior  accuracy  of  the  synergistic  model  as 
compared  to  the  life  prediction  using  linear  damage  summation. 

The  results  of  the  test  of  the  surface  flaw  specimen,  Figure  106,  indicate  that  the  total  life 
prediction  is  more  accurate  using  SINH  linear  damage  summation.  However,  the  accuracy  of 
the  surface  flaw  stress  intensity  factor  solution  diminishes  for  2b/w  >  0.45  due  to  the  prox¬ 
imity  of  the  edges  and  back  surface  of  the  specimen.  For  the  surface  flaw  specimens  tested  in 
this  program,  this  corresponds  to  a  half  crack  length  of  approximately  5.3  mm  (0.210  in.). 
Therefore,  the  data  of  crack  length  greater  than  this  value  are  not  entirely  valid  and  should 
not  be  used  for  model  evaluation.  These  data  are  shown  as  solid  triangles. 
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Figure  104.  Waspaloy  (PWA  1007)  Crack  Propagation,  Comparison  of  Data 
Generated  in  a  Surface  Flaw  Specimen  With  Baseline  Data, 

R  =  0.05,  0.167Hz  (10  cpm),  649°C  (1200° F) 
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Figure  105.  Waspaloy  (PWA  1007)  Model  Description  Test,  Mission  1, 
Compact  Specimen  No.  1565,  621°C  (1150° F) 
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Model  demonstration  tests  of  Waspaloy  fatigued  according  to  Mission  2  (Figure  101) 
were  also  predicted  well  by  the  SINH  model.  Test  results  from  the  compact  specimen.  Figure 
107,  revealed  good  accuracy  for  both  the  SINH  linear  damage  and  the  synergistic  SINH 
models.  Early  in  the  life  of  the  specimen  the  prediction  using  SINH  linear  damage  methods 
corresponds  closely  to  the  actual  fatigue  data,  while  the  synergistic  prediction  is  superior 
when  considering  the  total  life  of  the  specimen. 

Both  synergistic  and  SINH  linear  damage  calculations  of  crack  growth  resulting  from 
Mission  2  loading  were  less  adequate  in  predicting  the  life  of  the  surface  flaw  specimen. 
These  predictions  and  the  actual  data,  are  presented  in  Figure  108.  The  three,  inconsistently 
long,  crack  length  measurements  which  are  highlighted,  are  the  work  of  an  inexperienced 
technician.  Discounting  these  isolated  data  points  and  recalling  that  the  stress  intensity 
factor  solution  for  this  specimen  loses  accuracy  for  half  crack  length  greater  than  approxi¬ 
mately  5.3  mm  (0.210  in.),  both  life  calculations  are  demonstrated  as  conservative.  The 
approximate  values  Nca|c/Nac1ua|  for  the  linear  damage  and  synergistic  models  are  0.58  and 
0.73  respectively  for  the  5.3  mm  crack. 

b.  IN  100  (PWA  1073) 

The  IN100  demonstration  specimen  life  predictions  of  Mission  1  (Figure  101)  calculated 
using  the  interpolative  SINH  models  were  less  accurate  than  the  Waspaloy  predictions.  Both 
synergistic  and  SINH  linear  damage  calculations  gave  anticonservative  life  predictions  of 
both  specimen  geometries,  tested  at  710°C  (1310°F),  as  shown  in  Figures  109  and  110.  Test 
results  of  the  compact  specimen  (Figure  109)  show  the  synergistic  and  linear  damage  calcu¬ 
lations  to  overpredict  the  actual  specimen  life  as  demonstrated  by  the  respective  values  of 
Nc/Na  of  3.76  and  2.94.  The  prediction  errors  for  the  surface  flaw  specimens  were  similar 
(synergistic  model:  Nc/NA  =  1.57;  linear  damage  summation:  NC/NA  =  1.25). 

The  capability  of  the  IN100  interpolative  SINH  model  demonstrated  much  more  effective 
in  predicting  crack  growth  in  the  Mission  2  (Figure  101)  loading  sequence.  In  test  at  691  °C 
(1275°F)  of  a  compact  specimen  (Figure  1 11)  the  synergistic  model  prediction  was  extremely 
accurate  (Nc/NA  =  0.96),  while  the  life  calculation  using  methods  of  linear  damage  summa¬ 
tion  indicate  a  much  shorter  life  (NC/NA  =  0.64).  In  a  similar  Mission  2  model  demonstration 
test  using  a  surface  flaw  specimen  (Figure  112)  the  synergistic  model  prediction  exhibited 
reasonable  accuracy  (NC/NA  =  0.67),  while  the  linear  damage  summation  calculation  was 
significantly  conservative  (NC/NA  =  0.36). 
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Figure  107.  Waspaloy  (PWA  1007)  Model  Demonstration  Test,  Mission  2, 
Compact  Specimen  No.  1563,  621°C  (1150° F) 
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Figure  108.  Waspaloy  (P\V.  I  1007)  Model  Demonstration  Test,  Mission  2, 
Surface  Flaw  Specimen  So.  1498,  621°C  (1150°F) 
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Figure  110.  1N100  (PWA  1073)  Model  Demonstration  Test,  Mission  1, 
Surface  Flam  Specimen  No.  1473,  7 10°C  (1310° F) 


r 

141 


CROCK  LENGTH, 


TOTHL  MISSIONS 

PD  177371 

Figure  111.  INIOO  (PWA  1073)  Model  Demonstration  Test,  Mission  2, 

Compact  Specimen  No.  1333,  691°C  (1275° F) 
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Figure  112.  IN  100  (PWA  1073 )  Model  Demonstration  Test,  Mission  2, 
Surface  Specimen  No.  1575,  691°C  (1275° F) 
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S.  Critique 


The  basic  goal  of  the  model  demonstration  program  was  to  determine  the  effectiveness 
of  the  crack  propagation  models  to  predict  elevated  temperature  fatigue  crack  growth  in 
Waspaloy  (PWA  1007)  and  IN  100  (PWA  1073)  under  load  sequencing  representative  of  tur¬ 
bine  disk  operation.  This  evaluation  tests  the  capability  of  the  three-part  specimen  life  pre¬ 
diction  procedure  composed  of  the  segregation,  interpolation,  and  computation  algorithms. 

An  assessment  of  the  results  of  the  model  demonstration  program  may  be  accomplished 
by  constructing  a  log-normal  probability  plot  of  values  of  N  ,iCU|ated/NActuii|  for  the  collection  of 
demonstration  specimens.  Such  a  plot  is  presented  in  Figure  113.  Values  of  Nc/NA  which  are 
less  than  1.0  indicate  conservative  life  predictions  while  values  greater  than  1.0  are  anticon¬ 
servative,  the  statistical  analysis  of  the  test  results  provided  by  this  plot  shows  the  mean  of 
all  values  of  Nc/NA  to  be  1.07.  That  is,  the  average  life  prediction  given  by  the  synergistic 
model  is  7%  greater  than  the  actual  specimen  life.  This  highly  accurate  result  demonstrates 
the  excellent  capability  of  this  model  to  describe  cumulative  damage  effects  on  crack  growth 
in  Waspaloy  and  IN100.  This  may  be  contrasted  with  a  similar  analysis  of  results  of  the 
SINH  Linear  Damage  Model.  The  mean  value  of  these  predictive  errors  is  20%  less  than  the 
actual  specimen  life.  While  the  error  associated  with  the  nonsynergistic  model  was  not 
excessive  for  these  demonstration  tests,  a  considerably  more  significant  load  sequence  effect 
may  be  expected  under  more  representative  missions.  In  order  to  permit  usage  of  available 
testing  machines,  the  model  demonstration  missions  were  of  limited  cycle  count.  More  repre¬ 
sentative  mission  load  sequences  have  less  frequently  occurring  overloads,  resulting  in 
increased  crack  growth  synergism  (greater  error  in  the  nonsynergistic  life  calculation). 

While  the  average  model  behavior  is  represented  by  the  statistical  mean,  the  distribution 
of  the  sample  of  values  of  NC/NA  about  the  mean  provides  an  assessment  of  the  variability  in 
predictive  assuracy.  For  example,  the  plot  of  the  synergistic  model  results  (Figure  113)  indi¬ 
cates  that  for  one  hundred  future  life  calculations,  the  maximum  prediction  error  should  be 
anticonservative  by  a  factor  of  approximately  2-2.  One  in  1000  future  predictions  may  be 
anticonservative  by  a  factor  of  3.0.  Similar  results  occur  for  an  analysis  of  the  SINH  Linear 
Damage  Model  predictions. 

The  above  worst  case  analysis  considers  the  Waspaloy  and  IN100  data  together.  How¬ 
ever,  the  prediction  errors  for  Waspaloy  were  generally  less  than  for  IN100.  The  maximum 
and  minimum  errors  were  observed  in  IN  100.  The  primary  source  of  this  defference  is  not 
believed  to  be  the  difference  in  material.  Rather,  a  difference  in  testing  temperature  for  the 
two  materials  appears  to  be  responsible.  All  Waspaloy  demonstration  testing  was  conducted 
at  621°C  (1150°F),  while  the  Mission  1  and  Mission  2  IN100  tests  were  run  at  710°C  (1310°F) 
and  691°C  (1275°F)  respectively.  At  the  lower  temperture,  Waspaloy  displays  relatively 
limited  creep  crack  growth;  however  crack  propagation  in  IN100  under  sustained  load 
becomes  more  appreciable  in  the  indicated  temperture  ranges  for  these  tests.  Both  demon¬ 
stration  missions  have  significant  load  dwell  components,  and  at  the  IN100  operating 
temperatures  cracking  tends  to  violate  Linear  Elastic  Fracture  Mechanics.  Since  variability 
in  creep  crack  growth  data  is  generally  much  more  severe  than  observed  by  cyclic  fatigue 
crack  growth,  it  is  not  suprising  that  the  IN100  results  are  more  extreme  than  the  Waspaloy 
findings. 
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Figure  113,  Log-Normal  Probability  Plot  of  Synergistic  Model  Demonstra¬ 
tion  of  Test  Results 


In  order  to  examine  the  hypothesis  that  test-to-test  variability  in  crack  growth  was 
responsible  for  much  of  the  scatter  in  the  INlOO  model  demonstration  results,  a  number  of 
replicate  tests  were  performed.  The  most  severe  creep  conditions  were  selected:  Mission  1  and 
710°C  (1301°F).  The  results  of  five  such  tests  are  presented  in  Figure  114.  The  variability  of 
these  data  is  significantly  more  severe  than  observed  for  the  collection  of  the  separate  dem¬ 
onstration  tests  of  Figure  1 13.  The  extreme  data  scatter  for  the  replicate  tests  indicates  that 
much  of  the  variability  in  values  of  N*/NA  for  the  original  eight  specimens  was  the  result  of 
test-to-test  material  scatter  rather  than  an  arbitrary  error  in  the  cumulative  damage  model 
prediction. 

Another  source  of  model  prediction  error  was  extrapolation  of  AK  to  values  beyond  the 
limits  of  the  data  base  used  in  model  development.  The  range  in  AKover  which  data  ws 
generated  in  the  initial  test  program  wsa  necessarily  limited  by  the  scope  of  the  contract. 

While  it  is  difficult  to  assess  the  relative  contribution  of  errors  in  Cumulative  Damage 
Model  and  intrinsic  material  data  scatter,  the  combination  of  these  errors  is  represnted  by 
model  demonstration  results  of  Figure  1 13.  The  capability  of  the  synergistic  model  to  predict 
crack  propagation  under  loading  similar  to  demonstration  missions  1  and  2  is  described  by 
this  plot.  As  stated  earlier,  these  missions  were  chosen  to  be  extreme  cases  which  would 
exercise  the  model.  Typical  operating  missions  for  turbine  disks  are  expected  to  be  less 
severe,  therefore,  more  accurate  life  predictions  are  expected  for  field  usage  than  was 
observed  for  the  demonstration  test  program. 
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SECTION  III 
CONCLUSIONS 


The  primary  accomplishment  of  this  contract  was  the  development  of  a  computer  based 
mathematical  model  for  prediction  of  elevated  temperature  crack  propagation  of  Waspaloy 
(PWA  1007)  and  GATORIZED  *  IN  100  (PWA  1073)  subject  to  load  sequencing  representative  of 
turbine  disk  operation.  The  following  conclusions  are  drawn  from  the  results  of  testing  and 
analysis  which  supported  construction  of  the  model. 


1.  Loading  spectra  representative  of  turbine  disk  operation  contain  major 
and  minor  load  excursions,  as  well  as  load  dwells,  which  lead  to  synergis¬ 
tic  effects  on  crack  propagation.  In  the  missions  evaluated,  the  occurrence 
frequency  of  major  load  excursions  (overloads)  was  high  (generally  every 
50  cycles,  or  less),  and  the  observed  overload  ratio  (P()1/Pmai)  was  gener¬ 
ally  less  than  1.5.  Compressive  stress  excursions  may  occur  in  some  disks 
at  some  location  (eg.  cooling  holes). 

2.  Crack  propagation  data  (da/dN  vs  AK),  generated  under  constant  load 
amplitude  cycling  (and  repetitive  load  sequencing)  of  Waspaloy  (PWA 
1007)  and  IN100  (PWA  1073),  are  effectively  described  by  an  equation  of 
the  form 

fifi 

log  =  0.5  SINH  (C2  (log  AK  +  C.,))  +  C4 

where  C2,  C:J,  and  C4  are  empirical  functions  of  the  parameters  of  compo¬ 
nent  operation. 

3.  For  the  given  material,  Waspaloy  (PWA  1007)  or  IN100  (PWA  1073), 
equations  of  the  form  given  above  provide  the  capability  for  multiparame¬ 
ter  (eg.  frequency,  stress  ratio,  temperature,  overload  ratio)  interpolations 
to  define  representative  crack  growth  relationships. 

4.  Of  the  two  superalloys  Waspaloy  (PWA  1007)  and  IN100(PWA  1073),  the 
former  is  considerably  more  prone  to  test-to-test  variability  of  crack  prop¬ 
agation  properties.  The  basis  of  this  variability  is  the  large  range  of 
acceptable  crystalline  grain  size  for  W aspaloy  (PWA  1 007)  as  compared  to 
IN100  (PWA  1073). 

5.  U nder  constant  load  amplitude  cycling  with  dwell  at  peak  load,  the  rate  of 
crack  propagation  of  Waspaloy  (PWA  1007)  increases  with  increasing 
period  of  load  dwell.  However,  for  a  period  of  load  dwell  of  120  sec,  crack 
propagation  in  this  alloy  is  significantly  more  rapid  at  649°C  (1200°F) 
than  at  732°C  (1350°F).  The  higher  temperature  apparently  promotes 
creep-induced  crack  tip  blunting  which  reduces  the  rate  of  crack  growth  in 
this  alloy. 

6.  Negative  stress  ratio  fatigue  of  Waspaloy  (PWA  1007)  and  IN100  (PWA 
1073)  in  the  range  of -1.0  <  R  <  0.0  appears  to  produce  slightly  more  rapid 
crack  growth  than  observed  under  R  =  0  cycling.  This  effect  on  crack 
propagation  of  the  tensile-compreBsive  cycling  is  not  adequately  predicted 
by  an  extrapolation  of  the  effect  of  positive  stress  ratios. 
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7.  Elevated  temperature  crack  propagation  data  (da  dN  vs  AK)  of  Waspaloy 
(PWA  1007)  fatigued  at  representative  stresses  are  conservative  as  com¬ 
pared  to  similar  crack  growth  in  stress  fields  approaching  the  material 
yield  strength. 

8.  The  effect  of  prior  plastic  deformation  1.0%  strain  cycling;  tmMfl  =  on 

subsequent  elevated  temperature  (649°C)  crack  propagation  of  Waspaloy 
(PWA  1007)  is  to  increase  crack  propagation  slightly  above  characteristic 
cracking  behavior  of  virgin  (unstrained)  material.  In  IN  100  (PWA  1073) 
the  effect  of  similar  prestrain  is  neglible.  Test  results  of  both  materials  are 
within  the  material  scatterhand. 

9.  Under  repetitive  overload-underload-LCF  cycling  of  IN100  (PWA  1073)  at 
649°C  (12003F),  the  effect  of  the  compressive  load  cycle  following  the 
overload  is  to  reduce  the  beneficial  effect  of  the  tensile  overload. 

10.  A  periodic  dwell  (10  <  ANDwril  <  40)  at  peak  load,  interrupting  otherwise 
constant  load  amplitude  cycling  (R  =  0.1),  has  a  negligible  effect  on  crack 
propagation  of  Waspaloy  (PWA  1007)  and  IN100  (PWA  1073)  at  649°C 
(1200°F).  The  effect  of  similar  periodic  load  dwells  on  the  732°C  (1350°F) 
crack  propagation  of  Waspaloy  is  also  insignificant;  however,  the  periodic 
dwell  increases  crack  growth  rate  of  IN  100  at  this  temperature.  At  low 
values  of  AK,  this  increase  is  generally  predictable  using  linear  damage 
summation  methods. 

11.  Under  constant  load  amplitude  fatigue  (R  =  0.5)  interrupted  by  periodic 
overloads,  the  elevated  temperature  average  crack  propagation  rate  of 
Waspaloy  (PWA  1007)  and  IN100  (PWA  1073)  decreases  with  increasing 
overload  ratio  (1.0  <  OLR  <  1.5)  and  with  increasing  cycles  between 
overloads  5  <  ANol  <  40).  Over  this  range  of  representative  loading 
sequences,  the  effect  of  the  overloading  on  crack  growth  is  much  more 
pronounced  in  IN100  than  in  Waspaloy. 

12.  The  interpolative  model  of  elevated  temperature  crack  propagation  in 
Waspaloy  (PWA  1007)  and  IN  100  (PWA  1073)  are  effective  for  prediction  of 
cumulative  damage  effects  under  mission  loading. 
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APPENDIX  A 

OBSERVATION  OF  CRACK  RETARDATION  RESULTING  FROM  LOAD 
SEQUENCING  CHARACTERISTIC  OF  MILITARY  GAS  TURBINE  OPERATION 

J.  M.  Larsen,  C.  G.  Annls,  Jr. 


Pratt  &  Whitney  Aircraft  Group 
Government  Products  Division 
Mechanics  of  Materials  and  Structures 
P.O.  Box  2691 
West  Palm  Beach,  FL  33402 

ABSTRACT 


The  nature  of  crack  propagation  resulting  from  flight  loading  representative  of  military 
gas  turbine  operation  is  investigated.  Mission  stress  profiles  for  turbine  disks  fabricated  from 
the  superalloys  GATORIZED®  IN100  and  Waspaloy  contain  load  sequences  which  produce 
synergistic  effects  on  crack  propagation.  Major  load  throttle  excursions,  overloads,  occur 
routinely  during  flight,  and  a  retardation  in  subsequent  crack  propagation  generally  results. 
Such  mission  load  interaction  effects  have  been  addressed  in  crack  propagation  testing  employ¬ 
ing  repetitive  overload-fatigue  sequences.  The  influences  of  overload  ratio  (P  ,  ,/P  )  and 

the  number  of  fatigue  cycles  between  overloads  have  been  investigated  for  crack  propagation  at 
649°C  (1200°F),  and  an  interpolative  model  of  these  effects  is  presented.  A  determination  of  the 
instantaneous  crack  retardation  following  a  mission  major  load  excursion  is  accomplished 
with  an  unconventional  method.  The  existence  of  a  deceleration  in  crack  growth  rate,  delayed 
retardation,  following  a  mission  overload  is  verified.  Typically,  this  period  is  greater  than  the 
total  number  of  baseline  fatigue  cycles  applied  between  engine  mission  overloads,  and  delayed 
retardation  is  largely,  if  not  entirely,  responsible  for  the  beneficial  effects  of  the  overloading. 


Ksy  words:  crack  propagation,  fatigue,  retarding,  superalloyB 

Presented  at  ASTM  Symposium  on  “Effect  of  Load  Spectrum  Variables  on  Fatigue  Crack 
Initiation  and  Propagation”  San  Francisco,  C  A,  May  1979.  Submitted  to  ASTM  for  publication 
1979. 


INTRODUCTION 


Operational  loading  spectra  imposed  on  rotating  disks  in  a  military  gas  turbine  engine 
contain  load  sequences  which  differ  significantly  from  the  high  cycle  loading  encountered  in 
flight  of  an  airframe.  Fatigue  of  a  superalloy  engine  disk  is  a  low  cycle  phenomenon  resulting 
from  throttle  excursions  and  associated  thermal  stresses.  Major  and  minor  throttle  excursions 
compose  a  load  sequence  from  which  synergistic  crack  propagation  results,  and  this  load 
interaction  is  complicated  by  elevated  temperature  operation  and  concomitant  time  dependent 
behavior. 

Typical  loading  spectra  are  derived  from  service  missions  which  include  ferry,  training, 
and  terrain  following  radar  (TFR)  activity.  A  mission  composite  is  presented  in  Figure  A-l 
(Reference  1).  The  loading  spectra  may  conain  frequent  single  or  multiple  major  load  excursions 
(overloads)  with  a  small  number  (usually  no  more  than  50)  of  less  severe  throttle  excursions 
between  overloads.  As  a  result  of  this  frequent  overloading,  crack  growth  immediately  follow¬ 
ing  an  overload  is  of  increased  significance,  while  the  crack  retardation  commonly  observed 
many  cycles  after  an  overload  application  is  of  reduced  importance  since  this  retarded  growth  is 
soon  interrupted  by  another  overload  excursion. 

This  paper  reports  the  finding  of  an  investigation  into  the  effects  of  frequent  overloading  on 
elevated  temperature  fatigue  crack  propagation  (FCP)  of  the  super  alloys  G  ATORIZED  "  IN  100 
and  Waspaloy.  An  empirical  model  of  these  overload  effects  is  described,  and  the  characteristic 
retardation  of  crack  growth  which  follows  the  frequent  overloads  is  determined. 

CRACK  RETARDATION 

The  crack  propagation  resulting  under  the  complex  of  loading  and  environment  encoun¬ 
tered  in  a  military  gas  turbine  engine  is  not  predicted  satisfactorily  by  residual  stress  (Referen¬ 
ces  2  and  3)  or  closure  (Reference  4)  models,  and  a  need  for  further  study  of  this  isolated  problem 
exists.  The  current  work  has  employed  an  empirical  model  (References  5,  6  and  7)  of  elevated 
temperature  FCP  under  spectrum  loading.  The  basic  philosophy  of  this  approach  is  that  any 
complex  mission  spectrum  can  be  segregated  into  elemental  damage  events  which  can  be 
quantitatively  described.  The  crack  propagation  life  expected  under  such  a  spectrum  can  then 
be  computed  as  a  linear  addition  of  the  damage  associated  with  properly  segregated  events. 
Mission  segregation  is  based  on  the  definition  of  an  “elemental  damage  event”  as  the  smallest 
repeating  load-time  sequence  which  results  in  FCP  not  predictable  by  linear  damage  accumula¬ 
tion  alone.  A  simplified  mission  consisting  of  a  single  overload  followed  by  a  block  of  constant 
load  amplitude  fatigue  cycles  is  such  an  elemental  damage  event. 

Figure  A-2a  presents  a  schematic  of  a  crack  growth  curve,  crack  length  (a)  vs  number  of 
applied  fatigue  cycles  (N),  illustrating  the  common  effects  of  a  single  overload  on  previously 
unretarded  crack  growth.  When  the  crack  growth  produced  by  constant  load  amplitude  fatigue 
is  interrupted  by  an  overload,  crack  growth  accelerates  corresponding  to  the  overload  (Refer¬ 
ence  8).  Thereafter,  the  rate  of  crack  propagation  quickly  decelerates,  and  after  a  small  number 
of  subsequent  fatigue  cycles  the  rate  of  crack  growth  achieves  a  minimum.  This  deceleration  to 
a  minimum  rate  of  crack  propagation  is  known  as  delayed  retardation,  and  the  period  of 
delayed  retardation  is  defined  as  the  number  (NDR)  of  post-overload  fatigue  cycles  required  to 
achieve  the  minimum  crack  growth  rate.  Following  delayed  retardation,  crack  growth  con¬ 
tinues  at  a  near  minimum  rate  for  an  extended  period,  and  the  exhaustion  of  the  retardation 
process  is  marked  by  an  acceleration  of  crack  growth  to  regain  the  unretarded  rate.  The  total 
period  of  crack  retardation,  N*.  is  defined  as  the  number  of  cycles  during  which  crack  growth 
rate  is  retarded  following  an  overload  application.  The  nature  of  the  transient  crack  propaga¬ 
tion  behavior  which  results  following  an  overload  is  further  revealed  by  differentiating  the 
crack  growth  curve,  Figure  A-2a,  to  give  crack  propagation  rate,  da/dN,  as  a  function  of  N, 
Figure  A-2b. 
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Figure  A-2a.  Fatigue  Crack  Retardation  Resulting  from  the  A/ifilicath 


Figure  A-2b.  Fatigue  Crack  Retardation  Resulting  from  the  Application  of  a  Single  Oi'crloan 


Overloads  which  occur  at  cyclic  intervals  less  than  N*  interrupt  and  restart  the  crack 
retardation  process.  The  beneficial  effect  offered  by  intermittent  overloads  has  been  observed  to 
be  reduced  as  the  number  of  cycles  between  overloads,  DNOL,  is  reduced  to  fewer  than  N*  cycles 
(Reference  9  and  10).  As  D N0L  becomes  small,  the  relative  importance  of  delayed  retardation  is 
increased,  and  the  average  crack  growth  rate  associated  with  the  overload-fatigue  sequence 
increases. 


PROCEDURE 

Fatigue  crack  propagation  tests  of  the  nickel  base  superalloys  GATORIZED  IN100  and 
Waspaloy  were  conducted  on  electrohydraulic  testing  machines  operated  under  load  control. 
All  tests  were  performed  in  air  at  649°C  (1200°F)  and  the  loading  waveform  was  that  of  an 
isosceles  triangle  applied  at  a  frequency  of  0.167  Hz  (lOcpm).  Compact  specimens  (Figure  A-3) 
of  C-R  (circumferential-radial)  orientation  (Reference  11)  were  machined  from  disk  forgings 
representative  of  engine  disk  material.  The  test  specimens  were  pre-cracked  at  room  tempera¬ 
ture  using  standard  procedures  (Reference  11). 

Testing  was  designed  to  evaluate  the  influence  of  frequently  applied  periodic  overloads  on 
otherwise  constant  load  amplitude  ( AP)  fatigue  crack  growth.  A  schematic  of  the  basic  loading 
sequence  is  presented  in  Figure  A-4.  The  baseline  load  ratio  (R  =  minimum  load/maximum  load 
=  was  0.5  for  all  tests.  The  variables  of  test  were  the  overload  ratio  (OLR  =  (P0I/Pmu 

1.25, 1.50)  and  the  number  of  baseline  fatigue  cycles  between  overloads  (ANol  =  5,  20,  and  40). 

During  propagation  testing,  crack  lengths  were  measured  on  both  surfaces  of  the  specimen 
with  a  traveling  microscope.  To  facilitate  this  procedure,  the  test  was  interrupted,  and  the  mean 
load  was  applied  while  holding  the  specimen  at  the  temperature  of  testing.  The  increment  in 
crack  length  measurement  was  0.50  mm  (0.020  in.)  and  the  measurement  error  was  ±0.025  mm 
(0.001  in.).  The  resulting  "a”  vs  “N”  data  were  reduced  with  a  seven-point  incremental  polynom¬ 
ial  technique  (Reference  12)  to  produce  da/dN  vs  AK  data. 

RESULTS  AND  DISCUSSION 

The  crack  retardation  associated  with  a  single  periodic  overload  is  generally  reflected  in  a 
local  purturbation  in  the  “a”  versus  “N”  crack  growth  curve.  However,  as  ANol  becomes  small, 
the  resolution  in  crack  length  measurement  must  be  increased  in  order  to  perceive  the  transient 
variations  in  crack  growth  which  occur  between  successive  overloads.  In  the  present  experi¬ 
ment,  ANol  is  exceedingly  small  relative  to  the  increment  in  measurement  of  surface  crack 
length,  and  no  transient  variations  in  crack  growth  are  observed.  Rather,  the  collection  of  all 
“a”  versus  “N”  data  represents  a  trend  in  crack  propagation  resulting  from  very  frequent 
periodic  overloads. 

Reducing  the  “a”  versus  “N”  data  to  the  form  of  da/dN  versus  AK  data  and  performing  a 
least  square  regression  produces  analytical  functions  which  represent  the  average  crack 
growth  rate  for  the  specific  combination  of  OLR  and  DN0L.  The  regression  equation  was  based 
on  the  hyperbolic  sine  function  and  is  given  as:  (References  5,  6  and  7) 


log  (da/dN)  =  C,  sinh  (C2  (Jo(f  aR  +  ^  + 
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Figure  A -4.  Periodic  Overload  Fatigue 


The  increased  magnitude  of  the  overload  stress  intensity  range  (AKo1 )  is  not  reflected  in  the 
plotted  value  of  DK.  Rather,  the  overload  is  treated  as  an  isolated  variable  influencing  crack 
growth,  and  its  impact  is  revealed  in  the  dependent  variable,  da  dN.  Thus,  under  the  load 
sequencing  of  Figure  4,  DK  is  defined  as  K  -Kmin  for  all  cycles  (including  the  overload)  and 
da/dN  is  the  calculated  average  crack  growth  per  cycle  (including  the  overload). 

Defining  fatigue  with  periodic  overloading  as  a  succession  of  repetitive  overload-fatigue 
missions  (elemental  damage  events)  as  in  Figure  A-4,  one  may  plot  a  series  of  crack  growth  rate 
curves  for  the  collection  of  test  missions.  Figure  A-5  presents  such  data,  illustrating  the  effect  of 
aNoi  on  elevated  temperature  fatigue  crack  propagation  in  GATORIZED  IN  1 00.  For  purposes 
of  comparison,  a  single  regression  curve  (without  data)  representing  crack  growth  under 
constant  amplitude  loading  has  also  been  plotted.  This  baseline  crack  growth  is  observed  to  be 
essentially  equivalent  to  propagation  under  periodic  overloading  with  aN()|  =  5  (6  cycles  per 
overload-fatigue  mission).  The  equivalence  of  crack  growth  rates  results  from  a  balance  of 
overload  accelerated  crack  growth  and  the  subsequent  crack  retardation.  As  ANol  increases  to 
20  and  40,  the  average  effect  of  the  overloading  becomes  more  beneficial  than  damaging,  and 
the  resultant  propagation  rate  is  reduced  significantly  below  that  produced  under  constant  AP 
fatigue. 

The  dependence  of  da/dN  on  ANol  may  be  described  by  an  interpolative  model  of  the 
associated  SINH  curves.  The  coefficients  of  Equation  (Al)  are  defined  as  a  function  of  AN,,,: 

C{  =  a,  +  bt  (aN0L  +1)  i  =  2,  3,  4  (A2) 

Thus,  over  the  range  5  <  aN0L  <  40,  da/dN  is  given  by  a  continuous  function  (Equation  Al)  for 
which  the  three  coefficients  are  uniquely  defined  (Equation  A2).  Interpolations  on  ANol  define 
representative  SINH  curves  giving  average  crack  growth  associated  with  periodic  overload 
fatigue  of  ANq1  cycles  between  overloads.  This  permits  prediction  of  crack  growth  where  no 
actual  data  exist.  The  relationship  of  Equation  (A2)  also  provides  the  capability  for  limited 
extrapolation  beyond  the  data  base. 

The  effect  of  overload  ratio  on  fatigue  crack  propagation  may  also  be  described  by  an 
interpolative  model  of  the  SINH  coefficients.  Figure  A-6  illustrates  this  effect  for  periodic 
overloading  (ANol  =  40)  with  overload  ratios  of  1.0,  1.25,  and  1.5.  Data  of  overload  ratio  1.0  is 
generated  under  constant  AP  fatigue,  and  any  interpolative  model  of  OLR  must  converge  to  this 
condition  regardless  of  the  value  of  ANol.  The  expression  which  defines  the  SINH  coefficients 
as  a  function  of  overload  ratio  is 

Cj  -  dj  +  e,  (OLR)  i  =  2,  3,  4  (A3) 

Thus,  for  periodic  overload-fatigue  with  ANol  =  40,  the  combination  of  Equations  (Al)  and  (A3) 
uniquely  define  crack  propagation  for  1.0  <  OLR  <  1.5. 

Having  established  that  the  SINH  coefficients  are  linear  functions  of  overload  ratio  for  the 
case  of  ANol  =  40,  it  is  assumed  that  a  similar  linear  relationship  describes  the  dependence  of 
crack  propagation  for  ANol  <  40.  This  allows  combination  of  the  functions  (2)  and  (3),  permit¬ 
ting  full  interpolation  over  the  region  defined  by  1.0  <  OLR  S  1.5  and  5  S  ANol  <  40.  The 
coefficients  of  Equation  (Al)  are  given  by: 

C,  =  a(  +  fit  (OLR)  +  yt  (log  (ANol  +  1))  +  a,  (OLR)  (log  ( ANol  +  1 ))  (A4) 


i  =  2, 3, 4 
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Figure  A-5.  Effect  of  Cycles  Between  Overloads  on  Fatigue  of  IS  100,  Ol.R 
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Figure  A-6.  Effect  of  Overload  Ratio  on  Fatigue  of  IN  100,  AArn .  =  40,  0. 167  Hz,  649°  C 


This  Interpolative  SINH  Model  provides  capability  for  prediction  of  fatigue  crack  propagation 
in  IN100  as  a  function  of  overload  ratio  and  the  number  of  cycles  between  overloads.  This 
interpolative  model  is  limited  to  crack  retardation  effects  under  the  frequently  applied  over¬ 
loads  which  are  characteristic  of  military  gas  turbine  operation.  The  applicable  range  of  OLR 
and  ANol  may  be  extended  with  additional  testing. 

Development  of  a  similar  interpolative  model  of  crack  propagation  in  Waspaloy  has  also 
been  accomplished.  As  observed  in  Figure  A-7,  the  beneficial  influence  of  periodically  applied 
overloads  (OLR  =  1.5)  is  much  less  pronounced  in  this  material  than  in  IN100.  Crack  growth  in 
Waspaloy  with  ANol  =  5(6  cycles  per  overload-fatigue  mission)  is  considerably  more  severe  than 
propagation  under  constant  AP  fatigue,  and  the  average  effect  of  the  overloading  becomes 
beneficial  for  aNol  =  20  and  40.  In  spite  of  the  significantly  different  retardation  behavior 
between  Waspaloy  and  IN100,  an  interpolative  function  of  the  form  of  Equation  (A2)  is  satisfac¬ 
tory  representation  of  the  effect  of  ANol  on  crack  growth  in  this  alloy. 

The  model  of  the  effect  of  overload  ratio  (ANol  =  40)  on  crack  propagation  in  Waspaloy, 
Figure  A-8,  also  demonstrates  a  much  more  limited  effect  than  was  observed  in  IN100.  A 
relationship  of  the  form  of  Equation  (A3)  describes  SINH  coefficients  for  Waspaloy  crack 
propagation  as  a  function  of  overload  ratio.  Combining  the  models  of  the  effects  of  cycles 
between  overloads  and  overload  ratio  provides  interpolative  capability  for  both  variables 
simultaneously  (Equation  A4). 

The  coefficients  of  Equation  (A4)  are  given  in  Table  A-l  for  crack  propagation  of  IN  100  and 
Waspaloy  at  649°C  (1200°F),  0.167  Hz  (10  cpm),  and  R  =  0.5.  Substitution  into  Equation  (Al) 
gives  crack  propagation  rate  as  a  function  of  OLR  and  AN0L. 


POST-OVERLOAD  CRACK  RETARDATION 

An  Interpolative  SINH  Model  representing  crack  growth  under  frequently  applied  periodic 
overloads  is  an  effective  method  for  prediction  of  crack  growth  resulting  during  a  specific 
overload-fatigue  sequence.  However,  the  instantuneous  response  of  a  crack  to  the  overload- 
fatigue  block  is  not  evident.  From  knowledge  of  the  process  of  crack  retardation  resulting  from 
application  of  a  single  overload,  a  hypothesis  concerning  crack  growth  in  response  to  frequent 
periodic  overloading  may  be  established. 


The  crack  growth  curve  is  expected  to  range  between  two  general  extremes  depending  upon 
the  form  of  post-overload  crack  retardation.  For  a  given  value  of  ANol,  these  curves  are 
illustrated  in  Figure  A-9.  Assuming  delayed  retardation  does  not  occur,  the  minimum  crack 
growth  rate  should  immediately  follow  the  overload.  If  delayed  retardation  is  present,  the 
minimum  crack  growth  rate  should  be  observed  some  number  of  cycles  following  the  overload. 
If  ANol  is  less  than  the  period  of  delayed  retardation.  NDR,  a  subsequent  overload  interrupts  the 
delayed  retardation  process,  and  the  minimum  crack  growth  rate  should  immediately  precede 
the  overload. 
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Figure  A-9.  Fatigue  Crack  Propagation  Under  Periodic  Overload  Fatigue 
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Experimental  determination  of  the  nature  of  the  instantaneous  crack  growth  curve  in 
response  to  periodic  overloading  is,  however,  a  formidable  task.  As  previously  noted,  it  is 
extremely  difficult  to  obtain  accurate,  meaningful  measurement  of  crack  advance  between  very 
frequent  overloads  using  standard  optical  techniques.  Additionally,  measurement  of  crack 
length  at  the  specimen  surface  is  complicated  by  any  discontinuities  in  crack  advance  due  to 
microstructural  variations.  Such  discontinuous  fatigue  crack  growth  is  common  in  Waspalov 
and  to  a  lesser  extent  in  IN100.  Even  if  precise  measurement  of  the  surface  crack  length  were 
accomplished  without  microstructural  interference,  a  question  would  remain  concerning  the 
shape  of  the  crack  front.  That  is,  does  the  crack  front  geometry  remain  constant  during  periodic 
overloading,  or  do  periodic  changes  in  crack  curvature  accompany  the  overloads?  The  latter 
phenomenon  is  certain  to  complicate  data  interpretation.  Measurement  of  incremental  crack 
growth  by  determination  of  fatigue  striation  spacing  was  prevented  by  formation  of  a  heavy 
oxide  on  the  crack  face  during  elevated  temperature  testing. 

It  is  possible,  however,  to  determine  the  form  of  the  instantaneous  crack  propagation  curve 
indirectly.  Unlike  direct  measurement,  a  SINH  curve  describing  crack  propagation  under 
periodic  overloading  gives  average  behavior  over  the  period  of  many  overloads.  The  SINH 
representation  also  averages  crack  growth  discontinuities  due  to  variations  in  microsctructure, 
crack  front  geometry,  and  judgment  on  the  part  of  the  technician.  Thus,  each  curve  represents  a 
statistical  mean  of  the  microscopic  crack  growth  rate  as  AK  increases  with  crack  length. 

As  noted  earlier,  increasing  the  number  of  fatigue  cycles  between  successive  overloads 
produces  an  increase  in  crack  retardation  and  a  corresponding  reduction  in  average  crack 
growth  rate.  Alternately,  decreasing  AN()I  reduces  crack  retardation,  and,  in  the  limit  of  zero 
cycles  between  overloads,  the  da/dN  versus  AK  curve  should  approach  the  crack  growth  curve 
corresponding  to  constant  overload  fatigue.  For  the  current  example  of  constant  AP  fatigue 
(R  =  0.5)  interrupted  by  periodic  u”erloads  (OLR  =  Pu|  /Pmal  =  1.5),  the  parameters  defining 
constant  overload  fatigue  are  R  =  Pmm/P0|  =  0.33  and 

*^OL  “  ^bmelrne  [(OLR  -  R)/(l  -  R)J  -  SAK,,^  np  (A5) 

The  propagation  curve  for  R  =  0.333  fatigue  is  obtained  from  an  interpolative  model  of  the  effect 
of  load  ratio  (Reference  13).  Maintaining  the  convention  of  assigning  the  value  of  AKbaseline  to  all 
overload  cycles,  the  R  =  0.333  constant  AP  curve  is  translated  by  AKbMelinf  =  AKol/2  in  order  to 
represent  the  case  of  constant  overload  fatigue.  For  IN  1 00,  this  curve  is  presented  in  Figure  A- 1 0 
and  its  relationship  to  the  SINH  curves  representing  periodic  overload-fatigue  is  illustrated.  A 
similar  presentation  of  Waspaloy  data  appears  in  Figure  A-ll. 

From  the  collection  of  curves  for  a  single  material,  one  may  deduce  an  approximation  to  the 
instantaneous  form  of  post-overload  crack  propagation  due  to  periodic  overloading.  One 
assumption  is  required:  the  form  of  the  post-overload  crack  propagation  curve  is  similar  for 
ANol  =  5,  20,  and  40.  Noting  that  the  increase  in  crack  length  between  successive  overloads  is 
very  small,  the  associated  increase  in  AK  is  determined  to  be  correspondingly  small  and  of  little 
significance.  Calculating  crack  growth  rate  at  a  specific  value  of  aK  from  each  of  the  overload 
curves  gives  the  average  advance  per  cycle  under  the  indicated  loadir  ?  condition.  The  crack 
growth  during  each  of  the  individual  overload-fatigue  missions  is  then  calculated,  and  the 
difference  in  the  crack  growth  for  two  different  missions  defines  a  point  on  the  post-overload 
da/dN  vs  N  curve.  For  example,  by  subtracting  the  crack  advance  due  to  a  21-cycle  mission  from 
the  crack  advance  due  to  a  41-cycle  mission,  the  D a  corresponding  to  cycles  21  to  41  after  the 
overload  is  obtained.  The  average  value  of  da/dN  for  the  cycle  interval  of  21  to  41  cycles  after  the 
overload  may  now  be  obtained  by  dividing  Aa21  to4|  by  the  corresponding  AN  procedure  is 
given  by: 

da/dN21  to41  =  [41  <da/dN)4l  -  21  (da/dN)zl]/(41-21)  (A6) 
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Figure  A- 10.  Effect  of  Cycles  Between  Overloads  on  Fatigue  of  IX 100.  OI  K 


The  value  of  da/dN6lo2l  may  be  obtained  similarly,  and  da/dNlto6may  also  be  determined  by 
subtracting  the  standard  crack  growth  for  one  cycle  of  fatigue  at  the  overload  conditions  from 
the  6-cycle  mission  data. 

Plotting  the  calculated  values  of  da/dN  at  the  midpoints  of  the  cycle  intervals  (e.g.,  N  =  31 
for  the  interval  of  21  to  41  cycles)  gives  an  approximation  of  the  instantaneous  crack  propaga¬ 
tion  rate  following  an  overload.  For  IN  100  and  W aspaloy  at  AK  =  22  MPa  this  data  is  presented 
in  Figure  A- 12.  While  the  exact  path  of  the  immediate  post-overload  crack  propagation  rate 
curve  is  unclear,  as  indicated  by  the  broken  line,  the  existence  of  delayed  retardation  is 
confirmed.  In  each  alloy  the  crack  growth  rate  is  maximized  at  the  time  of  the  overload  and  the 
average  crack  growth  rate  decreases  in  each  of  the  subsequent  cycle  intervals,  i.e.,  1  to  6, 6  to  21, 
and  21  to  41. 


Figure  A-12.  Post  Overload  Crack  Growth  Exhibiting  Delayed  Retardation, 
OLR  =  1.5,  EK  =  22  MPa  /  m 


This  finding  is  verification  of  the  general  crack  behavior  which  was  hypothesized  earlier 
for  mission  crack  growth  subject  to  delayed  retardation.  That  is,  in  IN100  and  Waspaloy  a  crack 
subjected  to  constant  load  amplitude  fatigue  interrupted  by  periodic  overloading  will  display  a 
general  decrease  in  the  average  macroscopic  crack  propagation  rate,  but,  as  illustrated  in 
Figure  A-13,  the  instantaneous  response  of  a  crack  to  such  overloading  is  more  complex.  The  full 
extent  of  crack  retardation  does  not  develop  immediately  following  the  application  of  an 
overload;  rather,  post-overload  crack  growth  rate  decelerates  throughout  the  period  of  delayed 
retardation.  Typically,  this  period  is  greater  than  the  total  number  of  baseline  fatigue  cycles 
applied  between  mission  overload,  and  delayed  retardation  is  responsible  for  all  beneficial 
effects  of  the  overloading.  Since  crack  growth  rate  decelerates  during  the  entire  period  of 
delayed  retardation,  fatigue  cycles  immediately  following  an  overload  are  considerably  more 
damaging  than  later  cycles. 
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Macroscopic  Crack  Propagation 


FD  13233'A 

Figure  A-13.  Fatigue  Crack  Propagation  with  Frequent  Overloads 


This  behavior  deviates  significantly  from  the  general  long-term  effects  of  an  overload  on 
FCP;  therefore,  mission  crack  growth  in  IN100  and  Waspaloy  is  difficult  to  predict  by  any 
retardation  model  which  is  developed  solely  from  experimental  measurements  of  the  total 
period  of  retardation  N*.  Modeling  of  N*  addresses  fatigue  crack  retardation  as  a  steady-state 
process,  while  retardation  under  mission  loading  is  dominated  by  delayed  retardation  which  is 
a  transient  response.  Only  a  model  which  contains  this  transient  capability  is  expected  to 
effectively  predict  the  synergistic  effects  of  the  frequent  overloading  which  is  common  to  engine 
operation. 
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CONCLUSIONS 


1.  The  rate  of  elevated  temperature  FCP  of  IN100  and  Waspaloy  subjected  to  constant  load 
amplitude  interrupted  by  periodic  overloads,  decreases  with  increasing  overload  ratio  and 
number  of  cycles  between  overload. 

2.  The  combined  effects  of  overload  ratio  (1.0  <  OLR  5~  1.5)  and  number  of  cycles  between 
overloads  (5  <  ANol<  40)  on  crack  propagation  of  IN100  and  Waspaloy  are  effectively 
described  by  an  interpolative  model  based  on  the  hyperbolic  sine  function. 

3.  Post  overload  crack  retardation  resulting  from  frequently  applied  overloads  (OLR  =  1.5, 
AN0L<  40)  exhibits  delayed  retardation  in  IN  100  and  Waspaloy.  The  delayed  retardation 
is  not  exhausted  between  successive  overloads  and  is  entirely  responsible  for  any  reduc¬ 
tion  in  crack  growth  rate  in  these  alloys. 
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a  crack  length 

C;  SINH  material  coefficient  =  0.5 

C2  SINH  horizontal  scaling  coefficient 

C3  SINH  horizontal  inflection  coefficient 
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da/dN  cyclic  rate  of  crack  growth 
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Kmin  minimum  stress  intensity  factor 

Kol  overload  stress  intensity  factor 

N  number  of  cycles 

N*  number  of  cycles  during  which  crack  growth  rate  is  retarded  following  an  overload 

application 

Ndr  period  of  delayed  retardation 

P  applied  load 

Pmu  maximum  applied  load 

Pmm  minimum  applied  load 
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APPENDIX  B.  CUMULATIVE  DAMAGE  CRACK  PROPAGATION  TESTS 
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APPENDIX  C.  CRACK  PROPAGATION  DATA  FOR  SPECIMENS  TESTED  (Continued) 
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APPENDIX  C.  CRACK  PROPAGATION  DATA  I  OR  SPECIMENS  TESTED  (Continued) 
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APPENDIX  C.  CRACK  PROPAGATION  DATA  FOR  SPECIMENS  TESTED  (Continued) 
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APPENDIX  C.  CRACK  PROPAGATION  DATA  FOR  SPECIMENS  TESTED  (Continued) 
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APPENDIX  C.  CRACK  PROPAGA  TION  DATA  FOR  SPECIMENS  TES  TED  (Continued) 
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APPENDIX  C.  CRACK  PROPAGATION  DATA  FOR  SPECIMENS  TESTED  (Continued) 


1001335  PWA1073  1310F  MISSION  1001335  PUA1073  1310F  MISSION  1001473  PWA1074  1310  AIR  MISSION 

Rs  0-10  THICKNESSs12 .024  MM  R:  0.10  THICKNESS:12 . 824  MM  P  =  0.10  THICKNESS:  9.017  MM 

PMAX=  22.240  KN  HIDTH=63.269  MM  PMAX=  22.240  KN  UIDTH=63.269  MM  FMAX:  80.064  KN  WIDTH:23.660  MM 

COMPACT  SPECIMEN  COMPACT  SPECIMEN  SOLACE  FLAW  SPECIMEN 
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APPENDIX  C.  CRACK  PROPAGATION  I)ATA  KQR  SPECIMENS  TESTED  (Continued) 
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LIST  OF  SYMBOLS 


1  'nits 

Symbol 

Definition 

Metric  (Sit 

/'.ng/.u  h 

a 

crack  length 

mm 

in. 

B 

specimen  thickness 

mm 

in. 

ci 

SINII  material  coeffic  ient  =  0.’> 

dimensionless 

dimensionless 

Co 

SIXII  horizontal  scaling  coeflicient 

dimensionless 

dimen>i<  mless 

C3 

SIXII  horizontal  inflection  coeflicient 

log  ( AK  l 

log  ( AK ) 

C4 

SIXII  vertical  in  (lection  coeflicient 

log  (da/dX  1 

log  (da/dX) 

cc 

center  crack  specimen 

cs 

compact  specimen 

cl 

differential 

- 

da/dX 

cyclic  rate  of  crack  growth 

mm/cycle 

in. /cycle 

da/dt 

sustained  load  rate  ol  crack  growth 

mm /hr 

in. /hr 

A 

algebraic  difference 

Ad/ AX 

growth  rate  determined  from 

direct  secant 

mm/eycle 

in. /cycle 

AK 

applied  stress  intensity  = 

K  K 

rkmax  in 

MPas/  m 

ksiu  T7TT 

-A-'SOI, 

number  of  fatigue  cycles  applied 

between  periodic  overloads 

cycles 

cvcles 

^*^'d  well 

number  of  fatigue  cycles  applied 

between  period  load  dwells 

- 

((subscript ) 

final  value  of  parameter 

subscripted 

- 

f.g 

functions 

i(subscript) 

initial  value  of  parameter 

subscripted 

- 

K 

stress  intensity  factor 

MPasTirT 

ksi/  in. 

LCF 

low  cycle  fatigue 

- 

LKFM 

linear  elastic  fracture  mechanics 

- 

N 

number  of  cycles 

cycle 

cycle 

number  of  cycles  to  specimen 

failure  (actual) 

- 

Nc 

number  of  cycles  to  specimen 

failure  (calculated) 

- 

N* 

period  of  crack  retardation 

cycle 

cycle 

ndr 

period  of  delayed  retardation 

cycle 

cycle 

V 

test  frequency 

Hz 

Hz,  <  pm 

(cycles/minute) 

p 

load 

N 

Kips 

p 

max 

maximum  load 

N 

Kips 

Pmin 

minimum  load 

N 

Kips 

Pol 

magnitude  of  overload 

N 

Kips 

R 

load  ratio  =  P_  ■  /P_ 

min'  max 

- 

- 
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